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ABSTRACT 

We present the evolution of the volume averaged properties of the rest-frame optically luminous 
galaxy population to z ~ 3, determined from four disjoint deep fields with optical to near-infrared 
wavelength coverage. Our use of independent lines of sight substantially mitigates field-to-field varia- 
tions. We select galaxies above a fixed rest-frame y-band luminosity of 3 x 10^" h^Q Lq and charac- 
terize their rest-frame ultraviolet through optical properties via the mean spectral energy distribution 
(SED). To measure evolution we apply the selection criteria to a sample of galaxies from the Sloan 
Digital Sky Survey (SDSS) and COMBO-17 survey. The mean rest-frame 2200A through F-band SED 
becomes steadily bluer with increasing redshift but at all redshifts z < 3 the mean SED falls within 
the range defined by "normal" galaxies in the nearby Universe. Rest-frame 4000A/Balmer breaks are 
present in the volume averaged SED at all redshifts, indicating significant light and mass contributions 
from evolved populations, even at z ~ 2.8. We measure stellar mass-to-light ratios (A^*/L) by fitting 
models to the rest-frame UV-optical SEDs. Multiplying our volume averaged A^*/Ly estimates at 
each redshift by the measured mean V^-band luminosity density, we derive the stellar mass densities Pi,. 
The stellar mass density in galaxies selected at a fixed luminosity has increased by a factor 3.5 — 7.9 
from z = 3 to z = 0.1, where the range includes the uncertainty due to field-to-field variance within 
our own data. If we use our observed A^*/Lv evolution to select galaxies at a fixed mass, the stellar 
mass density evolves by a factor of 5.3 — 16.7. After correcting to total, the measured mass densities 
at z < 2 lie below the integral of the star formation rate (SFR) density as a function of redshift as 
derived from UV selected samples after a standard correction for extinction. This may indicate a sys- 
tematic error in the total or SFR(z) estimates. We find large discrepancies between recent model 
predictions for the evolution of p^, and our results, even when our observational selection is applied 
to the models. Finally we determine that Distant Red Galaxies (selected to have Js — Kg > 2.3) in 
our IV selected samples contribute 30% and 64% of the stellar mass budget at z ~ 2 and z ^ 2.8 
respectively. These galaxies are largely absent from UV surveys and this result highlights the need 
for mass selection of high redshift galaxies. 
Subject headings: Evolution — galaxies: formation 
content — galaxies: galaxies 
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Our knowledge of the high redshift Universe has in- 
creased rapidly over the past decade and is becoming 
ever more comprehensive. The initial large advances in 
this field were enabled by the efficient selection of z ^ 3 
starforming galaxies based on their optical light, e.g. 
iSteidcl ct al. (1996) and Stcidcl ct al. (1999), and by the 
development of efficient multi-object spectrographs on 
large 8-lOm class telescopes. Objects selected by these 
techniques allowed us for the first time to study statisti- 
cally significant samples of galaxies at z > 1.5. Resultant 
follow-up work on these optically (rest-frame ultraviolet; 
UV) selected objects at z ~ 3 (Lyman Break Galax- 
ies; LBGs) and subsequent samples at z ~ 2 (BM/BX 
objects; Steidel et al. 2004; Adelberger et al. 2004) 
have demonstrated that these galaxies have generally low 
extinctions E{B — V) ^ 0.15 and modest ages (< 1 
Gyr) and stellar masses {A4i,^ 10^°A^o) with a tail 
to hi gh end values exist ing at bright K magnitudes, 
e.g. iSawicki. Lin, fc Ye^ {1997), P apovich, Dic kinson, 
& Ferguson (2001; hereafter POD. i Shaplev et alJ l|200H) , 
.Shaple v et al. (2005), Rcdd v et all l)2005j) . 

Although efficient in telescope time for observing ob- 
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jects with R < 25, optical surveys will miss objects 
that have high extinctions or those with a lack of cur- 
rent vigorous star formation. In fact, most galaxies with 
A^*> 1O^^A^0 have optical magnitudes too faint for o p- 
tical spectroscopic foUowup ('van Dokkum ct al. 2006). 
A comprehensive view of the high redshift Universe re- 
quires that UV selection be complemented with selection 
from light at least as red as the NIR, which corresponds 
to the rest-frame optical out to z '--^ 3. Deep NIR obser- 
vations allow us to detect highly extincted (e.g. SCUBA 
galaxies; Small et al. 1997) or passive galaxies at z > 1.5 
(e.g. Daddi et al. 2005), which have very little rest-frame 
UV emission. 

However, getting deep NIR data is observationally ex- 
pensive and even over small fields requires substantial in- 
vestments of 8- 10m telescope time to reach depths fainter 
than Ks= 20 — 21. Until the recent development of 
meg apixel NIR a r rays p rogress was measurable but slow, 
e.g. IHogg et al.l l(l99l and IBershadv et alJ l|1998(l . A 
large step forward was taken with the Faint Infrared Ex- 
tragalactic Survey (FIRES; Franx et al. 2000; Rudnick 
et al. 2001; Labbe et al. 2003; Forster Schreiber et al. 
2005) and with the NICMOS observations of the Hubble 
Deep Field North (HDF-N; POl; Dickinson et al. 2003). 
FIRES is an ESO large program which imaged in Jg, H, 
and Kg the HDF-S and a field centered around the x-ray 
luminous cluster MS1054-03 at 2 = 0.83. At the time 
these fields had a unique combination of area and deep 
HST optical imaging and the ESO data extended this 
coverage to the i^T^-band. The FIRES data were particu- 
larly valuable because the associated deep if^-band data 
provide access to the rest-frame l/-band all the way out 
to z ~ 3. The very deep FIRES data allowed us to probe 
rest-frame U-band luminosities IV down to ~ 50% of the 
local Ly value. 

The powerful combination of very deep NIR and op- 
tical data now make it possible to assemble represen- 
tati ve samples o f the Univer se at z > 1.5. For exam- 
ple, iFranx et all (p003) and Ivan Dokkum et all 1)20031) 
have shown that galaxies selected to have (Js — Ks)> 
2.3 are luminous in the rest-frame optical, likely have 
high A^*, but would be almost completely missed by 
rest-frame UV selected surveys. Subsequent studies of 
these Distant Red Galaxies (DRGs) have shown that 
they have systematically older ages, higher extinctions, 
higher A^*, and comparable or even higher star forma- 
tion rates than the optically selecte d LBGs at compa- 
rable redshifts and K magnitudes ( van Dokkurn et aP 
200I iFdreter Schreiberltall I200I iLabbe et all 120051 



Papovich et alJl2005t Ivan Dokkiim et alJl200fih and may 

contribute significantly to the stellar mass density at 
z > 1.5 (Rudnick et al. 2003 - hereafter R03). Other 
author have also found galaxy populations that would be 
missed by rest-frame UV selection, highlighting the need 
for a comprehensive selection of galaxies (e.g. Daddi et 
al. 2004). 

While individual properties of galaxies are indeed im- 
portant, for some quantities the modeling of the vol- 
ume averaged population may yield more robust results. 
For example R03 demonstrated that bursty SFHs cause 
a smaller systematic error in the mean stellar mass-to- 
light ratio Aii,/L of the galaxy population - and hence 
estimates of Pi, - than when it is determined from av- 
erages of A^*/L estimates of individual galaxies. The 



utility of using cosmically averaged quantities has also 
been demonstrated by other authors across a range of 
topics. Using the luminosity density determinations at 
many wavelengths Madau. Pozzctti. & Dickinson (199^ 
tried to constrain the shape of the cosmic SFH by fitting 
the different bands simultaneously. Pei. Fall, fc Hausen 
([1999') derived independent constraints on the cosmic 
SFH by modeling the cosmic gas content, metallicity, 
and extragalactic background light. Using the FIRES 
data on the HDF-S R03 modeled the rest-frame optical 
mean colors and luminosity densities to derive the evolu- 
tion in the stellar ni ass d ensity at z < 3. At low re dshifts 
iBaldrv et all l|200l and iGlazebrook eTaLl l|200l mod- 
eled the mean galaxy spectrum from 2dF and the Sloan 
Digital Sky Survey (SDSS; York et al. 2000) to constrain 
the cosmic SFH. 

Using the deepest existing Kg-hand data R03 showed 
for the HDF-S that the evolution in the volume averaged 
rest-frame optical colors, e.g. {U — B)rest, {B — V)resti 
and {U — V)resti at z < 3 was monotonic, with bluer col- 
ors at higher redshifts. The colors also lay close to the 
locus of colors occ upied by individual normal galaxies in 
the local Universe ()Larson fc Tinslevlll978l: I.Tansen et alJ 
2000a). They also found that the evolution in the rest- 
frame colors could be approximated by smooth SFHs 
with a constant metallicity and dust. Models with these 
SFHs were then used to convert the [U — V)rest to 
7W*/Ly and hence to a mass density. 

A primary goal of this paper is to extend the analysis of 
R03. One major area of improvement is that this paper 
will model the full rest- frame UV through optical volume 
averaged SED of luminous galaxies in the Universe, al- 
lowing more freedom in the choice of dust attenuation, 
SFH, and age. The full rest-frame UV through optical 
information also gives us more insight as to the nature 
of the mean stellar population of luminous galaxies at 
high redshift. Perhaps the most important improvement 
with respect to R03 is that we include data from four 
disjoint fields, covering a total area of « 98.8 arcmin- 
utes. This is crucial since the field-to-field variations in 
the number densities of galaxies are expected to be large 
insmall fields. In addition, multiple fields are important 
to characterize the mean contribution of different galaxy 
populations, e.g. DRGs, to the volume averaged SED 
and Pi,. For example the field-to-field variations in the 
number densities of DRGs is known to be high with large 
differences between the HDF-N and HDF-S (Labbe et al. 
2003; hereafter L03). 

In this paper we combine FIRES data from the HDF-S 
and MS 1054-03 fields with those from the HDF-N catalog 
of POl and D03 and a catalog from the Great Observato- 
ries Origins Deep Survey (GOODS) images of the Chan- 
dra Deep Field South (CDF-S; Wuyts et al. in prepara- 
tion). Using photometric redshifts from these four cat- 
alogs we derive the rest-frame luminosities, luminosity 
densities and colors. We examine the trends in color and 
redshift and compute the full rest-frame 2200A through 
V^-band volume averaged SED. We then apply the same 
analysis to a galaxy sample from SDSS and from the 
COMBO-17 survey and use it to derive the evolution of 
mean SED properties in a consistent manner from z ^ 3 
down to = 0. Finally, we fit the full rest-frame mean 
SED with a set of simple models and use these to derive 
the mean yVf^^/Ly and p^, in each redshift bin. 
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In § 13 we describe the data. In § |3| we review the 
methods that we use to calculate photometric redshifts, 
rest-frame luminosities, luminosity densities, and global 
colors. In § ^ we present the results including the evo- 
lution of the rest-frame luminosity density, color, and 
volume averaged SED. In this section we also present 
the model fits to the mean SED and the derived A^^/L 
and Pi, evolution. We discuss these results in § |5l and 
present our conclusions and summary in § |3 Through- 
out this paper we assume VLyi = 0.3, JIa = 0.7, and Ho = 
70 hyo km s~^Mpc~^ unless explicitly stated otherwise. 

2. THE DATA 

We combine results from four different fields: the HDF- 
S, MS1054-03, HDF-N, and CDF-S. Below we briefly de- 
scribe the data reduction and the construction of cat- 
alogs. All magnitudes are quoted in the Vega system 
unless specifically noted otherwise. Our adopted conver- 
sions from Vega system to the AB system are Js,vega = 

Js,AB - 0.90, Hyeoa = Hab " 1-38, and Ks,vega = Ks_AB " 

1.86 ifBesseU fc Bre tt 1988). 

2.1. FIRES Data 

The reduced images and photometric catalogs for the 
two FIRES fields are available online through the FIRES 
homepage at 

http : //www . strw . leidenuniv . nl/ ~f ires. 

2.1.1. The HDF-S 

A complete description of the FIRES HDF-S obser- 
vations, reduction procedures, and the construction of 
photometric catalogs is presented in detail in L03; we 
outline the important steps below. 

The FIRES data on the WFPC2 field of the HDF- 
S is comprised of 101.5 hours of exposure with ISAAC 
(jlloorwood 1991} on the VLT. These exposures were 
split into 33.6, 32.3, and 35.6 hours in Js, H, and Ks 
respectively. The data were taken in service mode at 
the VLT and have a mean image quality better 0'.'48 in 
all bands. These NIR data were combined with the very 
deep o ptical data from WFPC2 taken by Cascrtano et al. 

Objects were detected in the Xs-ba nd image with ver- 
sion 2.2.2 of the SExtractor software l|Bertin fc ArnouO 
1199 61) . For consistent photometry between the space and 
ground-based data, all images were then convolved to 
0'.'48, corresponding to the effective resolution of the NIR 
band with the worst average seeing. Photometry was 
then performed in the [/300, B450, K06, -^814, Js, H, 
and ifj-band images using specially tailored isophotal 
apertures defined from the detection image. In addi- 
tion, a measurement of the total flux in the Ks-hand, 
Kl°^, was obtained using an aperture based on the SEx- 
tractor AUTO aperture, which includes a conservative 
aperture correction. The effective area of the HDF-S is 
4.74 square arcminutes, including only areas of the chip 
with exposure time > 20% of the total integration time 
in all bands. The accuracy of photometric redshifts (to 
be discussed in i i3.1(l and rest-frame optical luminosities 
in our ii's-selected sample is very dependent on the qual- 
ity of the NIR data (e.g. Rudnick et al. 2001) and once 
the S/N in the NIR bands becomes too low the quality 
of the photometric redshifts becomes too poor for a use- 
ful analysis. Since the Ks-ha.nd is our shallowest NIR 



band a cut there is a conservative proxy for a cut in the 
S/N of the other NIR bands, keeping in mind the range 
of galaxy NIR colors. The final catalog has 358 objects 
with ivr*°*< 23.14, which for point sources corresponds 
to a lOcr S/N in the custom isophotal aperture. 

2.1.2. The MS1054-03 Field 

The MS1054-03 observations, data reduction, 
and catalog construct i on ar e described in detail in 
iForster Schreiber et all l|2006D . This field has an x-ray 
detected cluster at z = 0.83 and at the time of the 
FIRES proposal was the field with the best combination 
of depth and area observed with WF PC2. The WFPC2 
data in Veoe and /814 are presented in lvan Dokkum et alJ 
(HOOO). The FIRES observations of this field consisted 
of 78 hours of Jg, H, and Ks imaging with ISAAC, 
supplemented with UBV imaging with FORSl on the 
VLT. The observations were split over four pointings. 
The total effective area is 23.5 square arcminutes, 
including only those pixels with exposure time > 20% 
of the total integration time in all bands. Observing 
conditions in the MS1054-03 field were generally similar 
to the HDF-S and since the exposures were split over 4 
pointings the depth is 0.7 magnitudes shallower than in 
the HDF-S. 

Object detection and catalog construction were per- 
formed in an identical way as in the HDF-S. The final 
catalog has 1380 sources with Kl°^< 22.34, which also 
corresponds roughly to a lOcr detection for point sources. 

2.2. HDF-N 

Th e WFPC2 data on the HDF-N are presented in 
Willi ams et all l|1996D . The data reduction of the NIC- 
MOS Jiio and iJiRn data and the Kitt Peak data 
i s des cribed by Dickinson (19!^ and [Dickinson et alJ 
(l200ri) as is the catalog construction^^. Objects were 
detected in a weighted sum of the Jno and -ffieo images, 
which are individually much deeper than the ii's-band 
image. The HST C/300, S450, Veoe, 4i4, Jiw, and Hieo 
data were then convolved to the same PSF and the fluxes 
were measured in matched apertures. The ground-based 
Ks data, which has much worse image quality, was ex- 
tracted using the "TFIT" method described in POl to 
achieve consistent colors between the space-based and 
ground-based images. A total magnitude in the Hiqq- 
band was estimated using the SExtractor AUTO aper- 
ture. The area of the HDF-N is 5 square arcminutes. 
The final catalog has 854 sources with H{fg{AB)< 26.5. 

The HDF-N is unique among our fields in that the the 
detection is not done in the if-band. Nonetheless, the 
i?i6o-band data arc sufficiently deep so that we are com- 
plete to the rest-frame V-band luminosity limit above 
which we select galaxies for this study (see § 13.2(1 . 

2.3. CDF-S 

From GOODS/EIS observations of the CDF-S (data 
release version 1.0) a fCs-band selected photometric cat- 
alog is constructed, described by Wuyts et al. (in prepa- 
ration). ISAAC imaging on the VLT provides 72 hours 
of exposure in J, 55 in H and 122 in Ks, split over 21, 

^•^ The reduced images and photometric catalogs are available 
from Mark Dickinson, mdickinsonOnoao . edu 
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12 and 23 pointings respectively. GOODS zeropoints 
were adopted for J and Kg. The 7J-band zeropoint was 
obtained by matching the stellar locus on a J — K vs 
J — H color-color diagram to the stellar locus in HDF-S 
and MS1054. The difference with the ofhcial GOODS 
77-band zeropoint varies across the CDF-S but on aver- 
age our iJ-band zeropoints is ~ 0.1 magnitudes brighter. 
ACS imaging provides photometry in i?45o, Veoe, -^814, 
and Z850 bands. All images were smoothed to match the 
NIR pointing with the worst seeing, 0'.'64. Identical pro- 
cedures as in the HDF-S and MS1054 fields were applied 
to detect objects and construct the catalog. A total ef- 
fective area of 65.6 square arcminutes is well exposed in 
all bands. The final catalog contains 1588 objects with 
Kl°^< 21.34 in this area. At Kl°*= 21.34 the median 
S/N in the Kg isophotal aperture is ~ 12. 

2.4. DRG selection 

Because the depth across our four fields is non-uniform 
we do not consistently apply a magnitude selection for 
DRGs. This may produce a slight error when compar- 
ing the DRGs to non-DRGs in each field. In addition, 
DRGs have different magnitude distributi ons in the dif- 
ferent fields ijForster Schreiber et al.ll2004|) and this may 
slightly bias the mass budgets determined in H4.4.1I In 
all fields our limiting Kg magnitude for DRG selection 
was dictated by the depth of our Js data, which must 
be deep enough so that non-detections in Js can still be 
selected as DRGs. 

In all fields DRGs were selected to have Zphot > 1-5. In 
the HDF-S 12 DRGs were selected to have Js-Ks> 2.3, 
Kl°^< 22.5 and Zphot > 1-5. In the MS1054-03 field 31 
DRGS were selected at i^*°*< 21.7. The brighter K, 
limit in comparison to the HDF-S is a direct result of 
the shallower data. In the CDF-S 65 DRGs were selected 
with Kl°'^< 21.14. This is 0.2 magnitudes shallower than 
the CDF-S catalog depth to insure reliable colors for ob- 
jects that are red in {Js —Kg). 

We select DRGs in the HDF-N to have Kf^K 21.05, 
which corresponds to a 5-ct detection in Ks- The very 
deep NICMOS Jno observations allow robust color mea- 
surements of red objects at this Ks limit. The selec- 
tion of DRGs in the HDF-N is complicated, however, by 
the use of the Jno filter, which is significantly bluer (ef- 
fective wavelength of l.lS/im for Jno and 1.25/im for 
Js) and about twice as broad (c.f. the NICMOS Instru- 
ment Handbook). It is therefore not appropriate to se- 
lect DRGs using a Jno — Ks > 2.3 color cut. Instead 
we compute a synthetic J, — Ks color using our best fit 
photometric redshift template (see ii3.1|) to derive the Jg 
magnitude. To assess the reliability of this method we 
first synthesize the Jno — Ks color and compare it the 
observed color. At Kl°^< 21.05 the outlier resistant nor- 
malized absolute median deviation {(Tnmad', equal to the 
rms for a Gaussian distribution) in the colors is 0.05 for 
all sources and 0.14 for those sources at 1.6 < z < 3.5. 
Computing the synthetic Js — Ks color for every galaxy 
we find two DRGs in the HDF-N with 1.6 < z < 3.5. 

3. METHODS 

3.1. Photometric Redshifts and Rest-Frame 
Luminosities 

Photometric redshifts Zphot for all galaxies are de- 
rived using an identical code to that presented in 



iRudnick et"all l)200lD and R03 but with a slightly modi- 
fied template set. This code models the observed SED us- 
ing non-negative linear combinations of a set of 8 galaxy 
templates. As in R03 we use the E, S b c, Scd , and Im 
SEDs from iColeman. Wu. fc Weed^mil] (11 98(1). the two 
least reddened starburst templates from iKinnev et all 
()1996ft and a synthetic template corresponding to a 10 
^yr year old simple stellar population (SSP) with a 
Saloctcr (1955) stellar initial mass function (IMF). In 
this paper we have added a 1 Gyr old SSP with a Salpeter 
IMF. 

Using spectroscopic redshifts Zspec we have determined 
that aNMAoizspec - Zphot / (1 + Zspec) = 0.06, 0.06, and 
0.08 at z < 4, < z < 1.5, and 1.5 < z < 4 respectively. 

It is important to remember, however, that the spec- 
troscopic samples in the HDF-S and HDF-N are highly 
biased toward UV bright objects at high redshift, e.g. 
LBGs. In the MS1054-03 and CDF-S field spectroscopy 
has been performed on Kg selected sources and those se- 
lected to be red in Js—Kg and in I—H. The photometric 
redshift accuracy for DRGs is cjnmad =0.11 and 0.09 in 
in MS1054 and the CDF-S respectively. There are, how- 
ever, very few DRGs with secure Zspec measurements (4 
in MS1054-03 and 3 in the CDF-S) causing the deter- 
mination of the Zphot accuracy itself to be uncertain at 
this point. From these first determinations, the accuracy 
of the Zphot measurements is still adequate for the con- 
struction of the luminosity densities and colors that we 
use here. 

The two sided 68% confidence intervals on Zphot, i-e. 
S Zphot , are computed using the monte carlo simulation 
described in R03. In general the differences between Zphot 
and Zgpec are well reflected by Szphot ■ 5zphot can therefore 
be used to judge the refiability of the Zphot estimates. 

Rest-frame luminosities are computed using the 
method described in the appendix of R03. We compute 
rest-frame luminosities 17^^^ in two UV bands centered 
approximately at 2200 and 2700A. These filters corre- 
spond to the ST - UV22 and ST - UV27 filters from 
(Bruzual & Chariot 2003; hereafter BC03). In addition 
we compute rest-fr ame luminosit ies in the UX^ B, V , 
R, and / filters of iBesselH l)1990j) . For these filters we 
use Mq^u = +5.66, Mq,b = +5.47, Mq.v = +4.82, 
Mq^u = +4.28, and Mq]u = +3.94. For' each galaxy 
we only compute U^^^ in rest-frame filters that are at 
shorter observed wavelengths than the Kg filter. In all 
cases where a spectroscopic redshift is available we com- 
pute the rest-frame luminosities fixed at Zgp^c- 

3.1.1. Star Identification 

Stars in all four fields were identified by spectroscopy, 
by fitting the object SEDs w ith stellar templates from 
iHauschildt. Allard. fc BaronI l|l999|) . and by inspecting 
their morphologies, as in R03. In well exposed regions 
of the HDF-S, MS1054-03, HDF-N, and CDF-S data 
we identified 29, 68, 16, 236 stars respectively down 
to Kl°*< 23.14, Kl°*< 22.34, Hlfa{AB)< 26.5, and 
if*°*< 21.34. All stars were excluded from the analy- 
sis 

3.2. Computing Luminosity Densities and Mean Colors 

In Figure we show the IV vs. cumulative enclosed 
co-moving volume and redshift for each field. The tracks 
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in each plot give the IV limit for a set of galaxy templates 
that corresponds to the observed 90% completeness limit 
in each image. From this plot we choose the minimum IV 
for which we are highly complete to the largest redshift 
over the largest area. The CDF-S image has the largest 
area but the data in this field are also our shallowest. 
Choosing an IV limit such that we are complete out to 
z < 3 in the CDF-S would cause us to throw away many 
galaxies with useful photometry in the other fields. In- 
stead we compromise by choosing a threshold such that 
we are complete out to z < 3 for the HDF-S, HDF-N, 
and MS1054, while still being complete out to z < 2.41 
in the CDF-S. This IV limit is L*^^'''''^ > 3 x 10^° h^^ Lq 
corresponding to the if*°*= 22.34 magnitude limit of 
the MS 1054-03 field. All computations of mean prop- 
erties presented hereafter are computed for galaxies with 
lV>Ly For the CDF-S, therefore, average quantities 
are only computed at z < 2.41. 

We compute the luminosity density i'^ every rest- 
frame filter in four bins which span the redshift ranges 
< z < 1.0, 1.0 < z < 1.6, 1.6 < z < 2.41, and 
2.41 < z < 3.2 and which are centered at redshifts 
of 0.73, 1.33, 2.01, and 2.8 respectively. These centers 
were chosen to split the co-moving volume in each bin 
evenly. The latter two bins have equal co-moving vol- 
umes. The 0.73 and 1.33 bins take this same volume and 
split it approximately 40/60%. This was done partly to 
sample better the large amount of time spanned by the 
z = — 1.6 interval and partly because of the rich clus- 
ter at z = 0.83 in the MS1054-03 field. We don't wish 
to bias our analysis by this rich cluster and by splitting 
the lowest redshift bin we exclude the lowest redshift bin 
(containing the cluster) in the MS1054-03 field without 
throwing out the valuable 1.0 < z < 1.6 information in 
that field. We compute j^^^^ in an identical fashion as de- 
scribed in R03, by adding up the luminosities of all galax- 
ies in each bin with lV>Ly"^°'*'^. A more sophisticated 
method such as a V/Vmax method is not needed since 
our threshold is in rest-frame luminosity and was chosen 
such that we are complete at all redshifts. As in R03 
we exclude galaxies with photometric redshift confidence 
intervals 5z/{l + z) > 0.16 but correct for their contri- 
bution to j™''* under the assumption that they have the 
same luminosity distribution as those sources with small 
redshift errors. Uncertainties in ji^'^''* are computed via a 
bootstrapping method in which 1000 samples are drawn 
from the subsample of galaxies with lV>Ly with re- 
placement allowed. The size of each bootstrap sample is 
a Poisson distributed number drawn from the measured 
number of galaxies with 1V>L^'^'"''^. For each iteration 
we recompute the correction for galaxies with large pho- 
tometric redshift confidence intervals and compute the 
luminosity density. We store the bootstrap iterations of 
the luminosity density in every band and use them later 
when performing the model fits in § 14.31 

It is important to note here that these j^^^^ values are 
lower limits to the total j^^^^, since we do not extrap- 
olate the luminosity function (LF) below Ly"'''\ We 

^■^ The MS1054-03 galaxies at z > 0.83 are lensed by the fore- 
ground cluster. We correct for the effects of lensing using the 
method described in Appendix El using the weak lensing map of 
this cluster by Hockstra ct al. (2000). Using the redshift distribu- 
tion of sources the inferred magnifications range from 5-20%. 



choose this method because the faint end slope of the 
rest-frame optical LF is observationally unconstrained 
at high redshift, for example, even our extremely deep 
i^s-band data in the HDF-S reach to only ~ 50% of the 
present day L*. To circumvent this limitation, many au- 
thors adopt a faint end slope from lower redshift bins 
and apply this to the higher redshift data. We prefer to 
avoid the uncertain extrapolation to the faint end of the 
LF by restricting ourselves to the observed data. 

To compute the luminosity densities averaged over the 
different fields we combine them weighting by their re- 
spective volumes, i.e. the total solid angle of each field. 
The mean colors we generate from the values, e.g. 

{U~V)rest = -2.5 X logioj[f^'+A/0,[; + 2.5 X logip jT'-^O.V- 

In such a formalism the mean color corresponds to the 
luminosity weighted mean colors of all galaxies with 
IV > L*^>'■'=''''. As in R03 we corrected the {U - B)rest 
for the {B — V)rest dependent contribution of emission 
lines, as determined using the spectra of the Nearby Field 
Galaxy Survey (NFGS; Jansen et al. 2000b). 

We point out that the measured values are weighted 
heavily to the largest field, which at z < 2.4 is the CDF- 
S. Because of the small number of fields and their very 
different weighting the determination of the field-to-field 
variance in our mean value of j"^* is not well statistically 
defined. 

3.2.1. Construction of Low Redshift Comparison 
Samples 

All deep surveys are limited by a small volume at 
z < 1, e.g. even the GOODS CDF-S field encloses only 
6.6 X lO^'Mpc^ or a cube of 40Mpc on a side, 50% of the 
distance to the Coma cluster. To build up a large volume 
at low redshift it is necessary to include a large redshift 
range in a single bin, which at low redshift smears to- 
gether large epochs of time. To supplement our data at 
these low redshifts we utilized data from the COMBO- 17 
survey and from SDSS. 

We computed at z < 0.2 for the SDSS Early Data 
Release (EDR; Stoughton et al. 2002) as in R03, de- 
riving j^"*'* in UBVRI from maps of number density as 
a function of L^"*^* and rest-fra me color. Th ese maps 
were computed as described in IBlanton et al . (2003a). 
While the Poisson errors in the SDSS are negligible, cos- 
mic variance and systematic errors do contribute to the 
uncertainties. For a more conservative error estimate, 
we adopt 10% errors on the luminosity density. For the 
SDSS luminosity function, our Ly '''''''^ encompasses 28% 
of the total light. 

To provide data at 0.2 < z < 0.9 we used 7^'^* measure - 
ments from the COMBO-17 survey ijWolf et all 12003(1 . 
which has a ^ 30 times larger survey area than the four 
deep fields combined. Specifically we used a catalog with 
redshifts of 29471 galaxies at z < 0.9, of which 7441 
had lv>Ly"°^'^ Using this catalog we calculated f^^^ 
identically to the deep fields. We divided the data into 
redshift bins of Az = 0.2 and counted the light from all 
galaxies contained within each bin which had lv>Ly 
The size of the redshift slices was chosen to well sample 

15 The J2003C catalog; available at 

http : //www. mpia.de/COMBQ/combo_index. html. 
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Fig. 1. — The distribution of rest-frame V luminosities as a 
function of enclosed co-moving volume and z is shown for MS1054- 
03, the HDF-S, the HDF-N, and the CDF-S in figures (a), (b), 
(c), and (d) respectively for galaxies with i^*°'< 22.34, Kl°'^< 
23.14, H*gjj(AB)< 26.5, and i<'*°*< 21.34. Galaxies which have 
spectroscopic redshifts are represented by solid points and for these 
objects L^"'^' is measured at Zspec- The large star in each panel 
indicate the value of L*^^^j from lBlanton et a l.' I'2003b') . The tracks 
represent the values of L^"''* for seven template spectra normalized 
at each redshift to the limiting magnitude in that field. The specific 
tracks correspond to the E (solid), Sbc (dot), Scd (short dash), Im 
(long dash), SBl (dot— short dash), SB2 (dot-long dash), and lOmy 
(dot) templates. The horizontal dotted line each panel indicates 
the luminosity threshold L^'''''''' above which we measure the rest- 
frame luminosity density j'^^^^ and the vertical dotted lines in each 
panel mark the redshift boundaries of the regions for which we 
measure jj^'^*' ■ 
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the 0.2 < z < 1 redshift interval while still being large 
enough to include many sources in each bin. The large 
solid angle of the COMBO-17 survey (0.75 deg^) assured 
that these slices still contained considerable co-moving 
volume. The formal 68% confidence limits were calcu- 
lated via a bootstrapping method as in ij l3.2l In addition, 
in Figure [21 we indicate the rms field-to-field variations 
a mong the three CO MBO-17 fields. As also pointed out 
in iWolf et aP l)2003j) , the field-to-field variations domi- 
nate the error in the COMBO-17 determinations. 

4. RESULTS 

In this section we first present our estimates of jy and 
of the evolution of the volume averaged color of luminous 
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galaxies. We then discuss the color evolution in terms 
of simple models. We also present the full rest-frame 
UV to optical volume averaged SED of luminous galaxies 
and show its evolution to higher redshifts, highlighting 
the mean SEDs of different subpopulations. Finally, we 
fit these mean SEDs with models and use the results 
to constrain the evolution in the global A^*/L and p^, 
also highlighting the contribution to the mass budget by 
different galaxy subpopulations. 

4.1. The Luminosity Density 

ROl were the first to estimate the evolution of the rest- 
frame optical luminosity function and luminosity density 
to z ~ 3, based on the number of rest-frame optically 
luminous galaxies in the HDF-S. In Figure |21 we show 
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and 



Ju as a function of cosmic epoch for our 



four deep fields and for the complementary z < 1 surveys. 
These estimates are given in Tableland Table As in 
R03, there is little evolution of jy^' at z < 3^^. FigureH 
also shows the large field-to-field variations present in 
our four deep fields. Large variations in jy among the 
three COMBO-17 fields are also seen, as indicated by 
the dotted error bars in Figure [3 Still, it is encouraging 
that the mean evolution in jy^^ is smoother than in the 
individual fields, indicating that we are getting closer to 
a measure of the true volume average value. 

The power of the COMBO-17 data is demonstrated 
here. With its inclusion it is apparent that there is a 
slight evolution toward brighter jy*** out to z ~ 1, as 
found by many authors (e.g. Lilly et al. 1996; Wolf et 
al. 2003). Moving to bluer wavelengths the evolution 
from z ~ 3 to z becomes more apparent, even within 
our own data. This already foreshadows the trends with 
color that we discuss in the following section. 

The field-to-field rms in jy^^ ranges from 20 — 75% 
in our 4 redshift bins. From this it is clear that data 
must be combined over large areas to robustly measure 
the evolution, especially for bright galaxies, which may 
be strongly clustered (e.g. Giavalisco & Dickinson 2001; 
Daddi et al. 2003). To mitigate field-to-field variations 
the most it is desirable for the individual fields to be 
larger than the typical clustering scale of objects in the 
field or to be spread over many independent lines of sight 
ifSomc rvillc et al. 2004). 

4.2. The Volume Averaged Color and SED 

We now proceed to measure the galaxy-averaged rest- 
frame colors at a given redshift/epoch. As discussed in 
R03 this average reflects a relatively smooth SFH for the 
ensemble of galaxies, even if individual SFHs are bursty. 
Therefore, a conversion of rest-frame color into a mean 
A^*/L will be more robust when modeling with a smooth 
SFH. In Figure |21 we illustrate this using an example 
which shows that the volume averaged color of a set of 
bursts behaves just like a population with the sum of 
their SFHs. 

In Figure 21 we present the redshift evolution of the 
luminosity weighted average color of all galaxies with 
IV > 3 X 10^° h^Q^ Lq, computed using Equation^ The 
colors become continuously redder from 3 to z ^ 0. 

^® The downward kink in j^^^^ at z ~ 2 is likely due to systematic 
effects in the photometric redshifts, which tend to preferentially 
depopulate the 1.6 < z < 2 region. 
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Fig. 2. — The evolution of 



J rest 



and j'u^* with time in 



panels (a), (b), and (c) respectively for galaxies with iV>L^''''''''. 
The solid symbols indicate the values for the individual fields as 
labeled while the large solid symbol is the mean value. The values 
are offset in time for viewing clarity. The errorbars reflect the 68% 
uncertanties in each field as characterized by a bootstrap simula- 
tion. The dotted errorbars represent the explicit rms field-to-field 
variations among our four fields. The open stars are the values for 
the COMBO-17 data and the open circle is the value for the SDSS. 
The dotted errorbars in the COMBO-17 data reflect the fleld-to- 
field rms uncertainties among the 3 COMBO-17 fields. Note the 
large field-to-field variation even between fields as large as 0.25 
square degrees. Note also that the MS1054-03 field was excluded 
from the lowest redshift bin because of the presence of a rich clus- 
ter and that we excluded the CDF-S from the highest redshift bin 
due to the depth of the NIR data in that field, which which made 
us highly incomplete at 1V>L^''''='' and z > 2.41. 
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Fig. 3. — A stylized example of how the volume averaged color of 
bursty SFHs reflects the sum of the individual SFHs. In the inset 
plot we show in the two lower panels two periodic SFHs with an 
identical duty cycle and peak SFR but which are perfectly out of 
phase. In the top panel we show the sum of their SFHs, which is 
equal to a constant. The main plot shows the evolution of the {U — 
V) colors associated with the individual SFHs and the evolution of 
a constant SFR model which is simply the sum of the individual 
SFHs (thick black line). We also show the {U — V) color which 
is derived from adding up the luminosities of each SFH in each 
band and using them to derive a total color as described in the 
text (stars). This luminosity weighted color is exactly the same as 
the color of a constant star formation model which is the sum of 
the individual SFHs. 



Contrary to the conclusions of R03, with our multi- 
ple fields it is apparent that the colors are only slightly 
less susceptible to field-to-field variations than the lumi- 
nosities, with the rms variations per redshift bin ranging 
from 5-50% for iU~V)rest and 20-75% for the This 
indicates that the mix of SED types also varies greatly 
from field to field, implying that large area surveys are 
even necessary to study the mean stellar populations of 
luminous galaxies in the Universe. 

In Figure we show our average color estimates and 
those from COMBO-17 and SDSS in color-color space. 
It is clear that all of the average estimates lie on a nar- 
row l ocus in this space. To in terpret this trend we re- 
fer to iLarson fc TinslevI l|1978j) who demonstrated that 
individual galaxies with normal morphological types lie 
on a narrow locus in the {U — B)rest vs. {B — V)rest 
diagram and that galaxies with large bursts of star for- 
mation scattered away from the relation, preferentially 
to blue ([/ — B)rest colors. To ascertain the amount of 
"burstiness" in our average colors we therefore compare 
them to those of individual local galaxies from the NFGS 
(Jansen et al. 2000a) in Figure |5l In general our col- 
ors and those of COMBO-17 are similar to the colors 
of local galaxies, implying a relatively smooth ensemble 
SFH. Our data in the highest redshift bin, however, have 
a slight systematic offset to bluer {U — B)rest colors at 
a given {B — V)rest color with respect to the local sam- 
ple. We explored whether part of this offset could result 
from systematic errors in the photometric redshifts. At 
z > 2 it appears that Zphot slightly overestimates the 
redshift on average, by about 0.2. These uncertainties 
likely result from an imperfect template set or an incor- 
rect amount of extinction. Although the redshift solution 
is driven primarily by the rough position of the domi- 
nant break in the SED, small color mismatches in the 
template set can slightly change the redshift. These sys- 
tematic redshift errors can indeed produce (U — B)rest 
colors which are bluer than the true values by up to 0.1 
magnitudes, although they can't entirely make up the 
differences with respect to the local data. The [U —B)rest 
color is especially sensitive to these effects, implying that 
more information than simply one color should be used 
to infer physical properties of the stellar population. The 
remaining differences between the high redshift and local 
samples can be understood, at least qualitatively, by con- 
sidering the effects of bursty SFHs on the mean color. Al- 
though the deviation from a smooth SFH is reduced when 
averaging over more galaxies there is still some residual 
burstiness left. 

We can also use the rest-frame UV to optical volume 
averaged SED to further constrain the nature of the mean 
stellar population. In Figure El we show the galaxy- 
averaged SED from the rest-frame UV through the opti- 
cal for our four redshift bins and compare it to the SEDs 
of local galaxy templates. The volume averaged SED at 
all redshifts falls well within the locus of normal local 
galaxies. Even in our highest redshift bin the mean SED 
is not as blue as a local starburst. Also evident is the 
presence of a break in the mean SED between the U and 
B bands. That this break is present even in our highest 
redshift bin indicates that the rest-frame optical light at 
z > 1.6 in rest-frame optically bright galaxies has signif- 
icant contributions from evolved stellar populations. 

In Figure we also split our sample into DRGs and 
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Fig. 4. — The redshift evolution of the volume averaged 
{U — B)reat and {B — V)reat- The symbols are identical to those 
in Figure|21with the addition that the {U — B)rest colors have also 
been corrected for probable emission line contamination as in R03. 
For clarity we have also plotted the errorbars of the individual 
fields as dotted lines. The colors of the lowest redshift COMBO-17 
point are artificially red due to the use of small central apertures to 
measure the color. For this reason we de-emphasize this point by 
only plotting a dotted errorbar. The integrated colors evolve red- 
ward with decreasing redshift with a much smoother progression 
than seen in the evolution of jy'^*' in Figure |^ The solid tracks in- 
dicate a set of solar metallicity exponentially declining SFHs with 
Zform = 4, E{B - V) = 0.25 (using a Calzetti et al. 2000 attenu- 
ation law), and with timescales of 1, 3, 6, 10, and 100 Gyr moving 
from top to bottom. Although these are merely example SFHs 
they highlight the inability of simple models with an exponential 
SFH, dust content, and metallicity to fit the time evolution of the 
global colors. 



Fig. 5.— The (U - B)reat vs. (B - V)re3t at z =0.73, 1.33, 
2.01, and 2.8 of all the stars in galaxies with 1V> 3 X lO^^ h:^^ Lq. 
The large symbols are identical to those in Figure]!] For clarity we 
have plotted the errorbars of the individual fields as dotted lines 
and we do not plot the ficld-to-field errorbars for the COMBO- 
17 data. The small solid points are the colors of nearby galaxies 
from the NFGS I Janscn et al. 2000ai) . which have been corrected 
for emission lines. The small crosses are the NFGS galaxies which 
harbor Active Galactic Nuclei. The thin tracks correspond to an 
exponentially declining SF H w i th a t imescale of 6 Gyr. The tracks 
were created using a ■ffalDeteil ^T95S^ IMF and the BC03 models. 
The dotted track has no extinction, the dashed track has been 
reddened by E{B — V) = 0.15, and the thin solid track has b een 
reddened by E{B - V) = 0.35, using the lCalzetti et alj |200C|) ex- 
tinction law. The thick black arrow indicates the reddening vector 
applied to the solid model track. The emission line corrected data 
lie very close the track defined by observations of local galaxies but 
are systematically bluer than all of the smooth models. 



non-DRGs and display their volume averaged SEDs. 
This is only done in the two highest redshift bins, cover- 
ing the redshift interval targeted, by design, by the DRG 
selection criterion. The DRG mean SED is uniformly 
redder than that of the non-DRGs. This is not surpris- 
ing at rest-frame optical wavelengths since these galax- 
ies were in fact selected to have red rest-frame optical 
(observed NIR) colors. The origin of the redder optical 
colors is partially elucidated by the detailed shape of the 
mean SED. It is apparent that the DRGs have stronger 
rest-frame optical breaks than the non-DRGs, indicative 
that evolved stellar populations are more prevalent in 
DRGs than in non-DRGs, to the same rest-frame opti- 
cal luminosity limit. Similar conclusions were reached by 
iForster Schreiber et al.l ("2004) and LOS based on the re- 
sults of SED fits to individual galaxies. DRGs, selected 
by their red rest-frame optical colors also have red UV 
SEDs. 

In summary, it appears in all cases that the volume av- 
eraged SED of all rest-frame optically luminous galaxies 
has a strong contribution from evolved stellar popula- 
tions and evolves to redder colors toward lower redshifts. 



4.3. Modeling the Volume Averaged Stellar Population 

In order to derive global Al^/L and subsequently p^, 
values from our mean SEDs it is necessary to use models 
to interpret them. In R03 we derived A^*/Ly estimates 
by using exponentially declining models applied to the 
global {U — V)rest color in which the SFR timescale, dust 
content, and metallicity were constant as a function of 
redshift With better data we attempt this same exer- 
cise again and find that the same simple models cannot 
accurately describe our higher precision data. We show 
fi ve examp le models in Figure ^ These models have 
a iSalpetej (llQSS*) IMF, solar mctaUicity, E{B - V) ^ 
0.25, Zform — 4 and exponentially declining SFHs with 
timescales of 1, 3, 6, 10, and 100 Gyr. Although these 
models represent only one possible set of parameters we 

In Figure 3 of R03 the tracks were mistakenly plotted with an 
extinction that varied with redshift. These tracks were correctly 
plotted with a constant E{B — V) in all other plots but the mistake 
in Figure 3 infiuenced our choice of best fit E(B — V) and therefore 
affected our absolute determinations of A4t /L and by about 0.1- 
0.2 dex. The main result of R03, the relative change in A^*/L and 
Pi, with redshift, however, was not affected since all values were 
computed with the same model. 
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Fig. 6.— The 2200A through /-band volume averaged SED of all 
the galaxies with iV > 3 X lO'^'^ \\^^ Lq. The different colored lines 
correspond to the different redshifts. The mean SEDs have been 
normalized to the same V-band rest-frame luminosity to highlight 
the evolution in the colors. The dotted lines are the SEDs of a set 
of galaxy templates, as labeled on the left side. The 2200A through 
7-band volume averaged SED at all redshifts falls well within the 
range of colors occupied by "normal" nearby galaxies. Even at the 
highest redshift bin the mean SED looks more like that of a late 
type galaxy than of a starburst. A break is also visible between the 
V and B bands, indicating the presence of an evolved population 
- with an associated 4000A/Balmer break - that is dominating the 
light. 



find it is generally impossible to simultaneously fit the 
([/ — B)rest and {B — V)rest colors, in the sense that 
the models have consistently too red {U — B)rest colors 
for a given (B — V)rest color for our data but with the 
opposite trend apparent for the COMBO- 17 data. This 
same disagreement is seen in Figure [S] where we plot our 
data, along with that of COMBO-17 and SDSS, in the 
{U — B)rest vs {B — V)rest color plane. The emission line 
corrections we applied to the colors did move the data 
in the direction of better agreement with the models but 
these corrections are already quite extreme at high red- 
shift (0.04 magnitudes in {U — B)rest) and larger shifts 
would be difficult to accommodate given the range in 
corrections inferred from the NFGS^^. To determine the 
robustness of the disagreement between the models and 
data we have explored different combinations of metallic- 
ity and extinction values and have found that no single 
set can reproduce these two colors at all redshifts. As we 
have mentioned already, however, the {U — B)rest color 
in particular is susceptible to redshift uncertainties and 
uncertainties in the emission line corrections. 
To mitigate the influence of one color on the global fit 

Of course, this is under the assumption that the EWs of emis- 
sion lines in the NFGS are the same as those of high redshift galax- 
ies at similar rest-frame colors. If, e.g. high redshift galaxies have 
higher equivalent widths then it would imply that we are under- 
correcting for emission lines. 



Fig. 7. — Same as for Figure ISl but splitting the two highest 
redshift bins split into DRGs (solid lines) and non-DRGs (dashed 
lines). At both redshifts the DRGs are redder than the non-DRGs, 
both in the rest-frame optical - where breaks are easily visible - 
and in the rest-frame UV. Because the DRGs still have light in the 
rest-frame UV the red colors are likely caused by a mix of old stars 
and dusty star formation. 



we choose to model the full UV-optical volume averaged 
SED to constrain the set of population parameters. As 
we have pointed out earlier, averaging over the whole 
galaxy population also averages over the SFHs of the in- 
dividual galaxies, making it more appropriate to apply 
smoothly varying SFHs when performing the model fits. 
This point of view is generally supported by the general 
agreement of our mean SEDs with those of "normal" 
local galaxies, which have been shown to be consistent 
with extended, relatively smooth SFHs (e.g. Kennicutt, 
Tamblyn, & Congdon 1994). With our choice of sim- 
ple models there are four possible free parameters that 
we consider: age, metallicity, dust, and SFH. We use the 
B C03 models an d an updated version of the Hyperz code 
jBolzonella et all i2b00|) to fit our data (Bolzonella pri- 
vate communication). We assume a Salpeter IMF with 
lower and upper mass cutoffs of 0.1 and 100 Mq respec- 
tively. Fitting with a different IMF, e.g. a Chabrier, 
would result in a scaling of the M^/L values by a con- 
stant factor at all redshifts (assuming that the IMF is 
universal). Thus the r elative trends would b e preserved. 
In all cases we use the iCalzetti et al.l ()2000j) attenuation 
law. We find that there exists some combination of pa- 
rameters that can fit the rest-frame UV-optical volume 
averaged SEDs as long as we don't require that the model 
parameters are the same at all redshifts. 

When fitting broad band photometry, especially only 
out to the rest-frame optical, there are very large de- 
generacies betwee n age, SFH, m e tallici ty, and dust con- 
tent, e.g. POl, iShanlev et alJ (2001). For example, 
in some cases high redshift SEDs can be fit reason- 
ably well with both a constant SFH (hereafter CSF), 
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moderate ages, and high extinctions and SFRs, or with 
an exponentially declining model with lower extinc- 
tion and SFRs. Despite these degeneracies, however, 
A^*/L are better constrained, since different combina- 
tions of stellar population parameters can give simi- 
lar d e rived mass-to-light ratios iForster Schreiber et al. 
120041: Ivan Dokkum et all l2004|) . We will therefore use 
our fits to measure but will not discuss the con- 

straints on the stellar population parameters in detail. 
There is tentative evidence from NIR spectroscopy of 
bright galaxies at z < 3 th at met allicitics ar e typically 
around solar i jvgii Dokku m et al.l [2004; Sh aplev et al.l 
I2004t Ide Mello etal.. .2004) and we therefore limit our 
fits to solar metallicity. 

To span the range of possible models we perform fits 
with either a CSF, Simple Stellar Population (SSP), or 
exponentially declining model with a timescale of 300 
Myr (rsooMyr)- In all cases we restrict the age of the 
population in each redshift bin to be younger than the 
age of the Universe for the middle redshift in that bin. 
We parameterize the amount of extinction by the mag- 
nitudes of attenuation in the l^-band Ay- For the CSF 
and TaooMyr models we allow the extinction to vary be- 
tween < Av < 3.0 and for the SSP models between 
< Ay < 0.5. This latter choice is to prevent the fitting 
from choosing highly extincted very young models with 
no ongoing star formation, which would result in unphys- 
ically large intrinsic luminosities. The upper limit of 3.0 
magnitudes of extinction does not affect our conclusions. 
For each redshift bin we choose the best-fit SFH, age, 
and extinction combination and use the AA^/hy of that 
model, which includes the stellar mass loss from evolved 
stars. We then multiply that A^^^/Ly by jy^*^ to derive 

To derive uncertainties on A^^^/Ly and we fit the 
ensemble of mean SEDs produced by our bootstrap sim- 
ulation described in § 13.21 Each bootstrapped SED was 
computed from the same set of galaxies, which is neces- 
sary because each galaxy contributes luminosity in every 
passband, i.e. the values are correlated. We then 
use the distribution of Al^/Lv^ and values from the 
fit to the ensemble of SEDs and determine the 68% con- 
fidence limits. 

The best fit SFH for each volume averaged SED is 
presented in Figure |H| We also indicate the formal best 
fit SFH, which is the raooMyr model for the three lowest 
redshift bins and the SSP model for the highest redshift 
bin. It is important to note, however, that the raooMyr 
model is statistically allowed at all redshifts and that our 
A^*/L determinations are not very sensitive to the exact 
SFH. As was visible in Figure El the rest-frame optical 
break is present at all redshifts, with the strength of the 
break decreasing toward higher redshifts. This is also 
reflected in the model fits. 

We also fit the mean SED from the SDSS survey to con- 
strain the evolution of Al^/L and to low redshift. The 
COMBO-17 data were very useful for characterizing the 
overall trends of color and for assessing the applicability 
of certain models for the color evolution. Nonetheless, 
we do not model the COMBO-17 data in detail since 
the lack of deep NIR observations limits us to only a 
relatively small range in rest-frame wavelength at each 
redshift, with correspondingly poor constraints on the 
best-fit Ai-t/Lv- The low redshift sample from the SDSS 



is crucial as our IV cut complicates the use of literature 
values for p^,, which are always quoted as total values. 
By fitting the mean SED of the local sample as deter- 
mined with our IV threshold we can consistently track 
the evolution over redshift. The only differences in ana- 
lyzing the low and high redshift samples is in the allowed 
stellar population parameters and the derivation of un- 
certainties. Because luminous galaxies from the SDSS 
are alm ost entirely made uP of evolved early type p opu- 
lations l|Blanton et alJl2003at Ikauffmann et alJl2003|l we 
limited the attenuation to 1 magnitude in the V-band. 
In addition it is difficult to obtain a realistic uncertainty 
estimate for the SDSS. The uncertainties in j^^^^, and 
hence in A4i,/h and are dominated by systematic un- 
certainties. Therefore we give a 10% error to all derived 
SDSS quantities. 

4.4. Ai^/Lv and 

In Figure El we show the evolution of the derived 
(A^*/Ly) for luminous galaxies as a function of time 
and redshift. These are also given in Table 01 As ex- 
pected the bluer colors of the volume averaged SEDs and 
the decrease in the break strength toward higher redshift 
implies a lower value of (A4+/Ly). (A4*/Ly) increases 
by a factor of ^ 10.7 from z = 2.8 to z = 0.1. 

Using the (M^/Ly) values and the jy measurements 
we derive p^ for galaxies at high redshift and in the 
SDSS with lV>Ly"'°^'^. Our absolute p^ measurements 
at < z < 3.2 are given in Table |31 The mass density in 
luminous galaxies between decreases by a factor of « 4.8 
between z — OA and z — 2.8 and ^ 50% of the stellar 
mass was in place by z > 1. This is in broad agreement 
with the results of R03, who measured a decline by a fac- 
tor of ^ 10 in p^, between z = 0.1 and z = 2.8, although 
with significantly larger errorbars and with no account- 
ing for cosmic variance. The rms field-to-field variations 
in p^ among our four fields ranges from 40 — 61%. 

To compare our data to that from other authors we 
scale our data to "total" values. We calculate an upwards 
correction factor of 2.1 by measuring the total SDSS p* 
value and comparing it to p^ for the SDSS galaxies with 
lV>Ly"'''^. After correcting all our measurements to to- 
tal we normalize our values to the z = p* measurement 
of COl and plot them along with those of other authors 
in Figure ^1 The relative values are given in Table |31 
We also determine how the inferred evolution in the total 
p^, changes if we allow Ly '''^''^ to evolve in the same way 
as our Mi^/Jjy values. The change this can induce is in- 
dicated by the vertical extent of the boxes in Figure [TUl 
and is reflected in the confldence interval in Ap^ given 
in Table El Our local total mass density measurement 
from SDSS is « 25% lower than that of COl. Some of 
this discrepancy may be due to the slightly higher mean 
redshift of our data. COl estimate their stellar masses 
at z = whereas the mean redshift fo r the SDSS sample 
is 0.1. Using the z = 0.1 SFRD from iBrinchmann et alJ 
ll200l we would predict that p* would increase by w 9% 
between z = 0.1 and 0, accounting for some of the dis- 
crepancy. Another possible source of the discrepancy is 
the use by COl of a fixed formation redshift of z = 20 for 
their SED fittin g derived A ij, estimates. 

According to IB ell et al] f2003) the old age assumed 
for all galaxies could bias their estimates by w 10% 
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Fig. 8. — The best fit models to the volume averaged SED. Each row corresponds to one of the redshift bins and each column corresponds 
a different SFH. The panel in each row with a star indicates the formal best fit SFH, although the other SFHs are generally allowed within 
the uncertainties. The filled points show the measured values and the open circles show the best fit model fluxes. As in Figure 151 the 
rest-frame optical break of both the data and best-flt model is apparent at all redshifts and increases toward lower redshifts. 



iBell et al.' (200^'^^. By refitting the SEDs of galaxies 
with SDSS an d 2MA SS data and letting the age vary 
iFbntana et al.l lj2004j) estimated that forcing an old age 
can bias the masses from COl upwards by < 20%. Tak- 
ing all these effects into account, the remaining differ- 
ence between our SDSS point and that of COl falls well 
within the range of expected systematic errors in our p^, 
estimates. 

Compared to the COl point, the mass density at z = 
2.8 is 5-8 times less than at z = 0, where the range cor- 
responds to our formal 68% confidence limits for a fixed 
j^thicsh^ As discussed in ti4.5l however, it may not be 
appropriate to compare the evolution of our IV limited 
sample to that of the total estimates presented by dif- 
ferent authors (see below). If we instead allow Ly"^"^*^ to 
vary along with our A^^/Ly values the mass density at 
z — 2.8 is constrained to be 5.3-16.7 times less than at 
z = 0. 

19 The determination of IBell et aF f200?) agrees excellently 
with that of COl and also assumes a fixed old age (12 Gyr) when 
deriving A4i, of galaxies in their sample. 



At higher redshifts is estimated by other authors 
using SED fits of individual galaxies. In all cases 
we have converted to a Sa lpeter IMF when necessary. 
iBrinchmann fc Ellis! l)2000|) fit a set of /-band selected 
galaxies at z < 1 from the Canada France Redshift 
Survey (CFRS). They are only complete over the range 
10.5 < logA^* < 11.6 and their value here includes a 
20% correction to a total mass d ensity re commended by 
IBrinchmann fc Ellis ( 2000). Drorv et all (^004) estimate 
Mi, for z < 1.2 galaxies with K < 19 over 0.28 deg^ from 
the Munich Infrared Cluster Survey (MUNICS). They 
also adopt an SED fitting technique but use their best- 
fit models to derive A4i,/hK and then extrapolate their 
SEDs from the observed to the rest-frame /C-band to de- 
rive IK. Using the DEEP 2 s pectroscopic su rvey coupled 
with K-hand imaging, Bundv et al.l l)2006j) use an SED 
fitting technique to estimate galaxy masses and then fit 
a Schechter function and extrapolate it to get their p^, 
estimates. At high redshifts/faint magnitudes, they sup- 
plemented their spectroscopy with photometric redshifts 
but these have a minimal effect on the derived values. 
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The errorbars for the lBundv et al.l l)2006j) po ints are from 
the va riance in pi, between their four fields. iBorch et al.1 
l)2006|) computed the stellar mass density from 25,000 
galaxies in the COMBO- 17 survey. They fit their combi- 
nation of 17 broad and medium optical bandpasses with 
SEDs to derive the stellar masses for each galaxy. They 
then fit Schechter functions to thei r data in each redshift 
bin and integrate this to derive p^. iFontana et al.l 1)20041) 
use the K20 survey to measure the stellar mass function 
out to z = 2 using stellar masses derived by SED fitting 
of individual galaxies. They then integrate this mass 
function to obtain p-^. In their 1.5 < z < 2 bin their 
mass function is complete only down to the ~ M.* value 
of their best-Schechter fit to the mass function and the 
extrapolation to lower M.i, is large and uncertain. D03 
and F03 derive at z ~ 3 from the HDF-S and HDF- 
N WFPC2 fields respectively, where the raw data in the 
HDF-S are in common with the FIRES data and where 
the HDF-N catalog is identical to what is used in this 
paper. The techniques in both these papers are nearly 
identical, using SED fitting of individual galaxies to de- 
termine the mean A^^/L for all galaxies and then apply- 
ing this to the luminosity density as determined by inte- 
grating a Sc hechter ftinction fit to the rest-frame optical 
luminosities. ' Drorv et al.l ((2005*) estimate the mass den- 
sity in the FORS Deep Field (FDF) using an /-band se- 
lection sample and in the GOODS-S field using a K-hand 
selected sample drawn from the same raw data as used 
for this paper. They derive mass densities by fitting the 
observed SEDs and then add up all the galaxies in their 
sample at each redshift. Although they measure p^, out to 
z ^ 5 we only plot it out to z 2 for two reasons. First, 
their /-band selection in the FDF corresponds to a selec- 
tion in the UV at z < 1.5 and the authors are therefore 
more sensitive to unobscured star formation than to stel- 
lar mass. They point out that their /-band data is deep 
enough to detect all but 10% of the /C-selected objects 
in the deep FIRES images of the HDF-S but it is exactly 
these missed objects which comprise a significant fraction 
of the mass density at high redshift, i.e. DRGs. Second, 
although they do have A'-band data in the GOODS-S 
field it is not deep enough to be mass complete for any 
but the most massive objects at z > 2, but they add up 
all objects, regardless of whether they are complete for 
those masses. It is therefore difficult to interpret their to- 
tal values at z > 2. Glazcbrook et al. (2004) determine 
stellar masses by fitting the observed SEDs with stellar 
population synthesis models and then compute using 
a V/Vmax method. The compute down to different 
mass limits and note their decreasing incompleteness in 
mass with increasing redshift and decreasing mass. After 
converting to a Salpeter IMF we plot their points down 
to their lowest mass limit, 10^° ""^ Al©. At the highest 
redshift bins they are incomplete and the resultant p* 
values can only be considered as lower limits. 

In our < z < 1.0 redshift bin our data fall slightly 
lower than the other determinations. This difference may 
be because the CDF-S is underdense at low redshifts. 
We suspect this because IWolf et all ()2003|) note that the 
Extended CDF-S is underdense at 0.2 < z < 0.4 and 
that this trend apparently is true at z < 0.6 (Papovich 
private communication) . Although there is a loose struc- 
ture in th e Extended CDF- S at z = 0.67 and a cluster at 
z = 0.73 lGini et alJ l|2003D . these should not dominate 



the counts over the much larger CDF-S field. We don't 
include the MS1054 data in our < z < 1 measurement 
because of the massive cluster in that field. However, 
the large field-to-field variation among our 4 fields at 
z < 1 makes our data consistent with those from other 
surveys. At z > 1 our p^, estimates agree very well with 
those of the other authors but with smaller formal er- 
rors corresponding to the lower uncertainties afforded by 
the multiple fields. The only notable discrepancy is with 
the HDF-N determinations of D03; we address a possible 
cause for this in the following subsection. 

We also show the curves that correspond to the in- 
tegral of the parametric fit to the SFR(z) curve from 
(Cole et al. 2001; COl), both with and without a sub- 
stantial extinction corrcction^°. As pointed out by, e.g. 
D03, R03, F03, and Fontana et al. (,2004), an extinction 
correction to the UV-dcrivcd star formation rates is re- 
quired to match the p^ measurements at all redshifts, 
although the nature of this correction is highly uncer- 
tain. The p^ evolution predicted from the UV selected 
extinction corrected SFR is in broad agreement with the 
direct p* measurements, notably at z ~ 3. Nonetheless, 
the measurements, including those from the literature 
are systematically below the extinction corrected curve 
by about 0.2-0.4 dex at z < 2, with the exception of the 
iDro rv et al. (2004) data. This slight offset could be due 
to an erroneously high extinction correction to the SFR 
estimates, but these can't be changed too much without 
coming into disagreement with the local p^, determina- 
tions. Additionally, some of the disagreement may come 
from systematic underestimates of p^, although it is not 
immediately clear what effects would plague all surveys, 
which use different techniques to estimate the mass den- 
sities. Finally, it may be that the galaxies which enter 
into the SFR(z) determinations are too faint in the rest- 
frame optical to enter the NIR selected samples. The star 
formation rate density (SFRD) measurements rely on an 
extrapolation to the faint end of the rest-frame UV lu- 
minosity function. Estimates of the faint end slope have 
a large range, with the original determination from the 
HDF-N by Steidcl ct al. (1999) being « -1.6 but with a 
later measureme nt from the FORS D eep Field giving a 
value of ~ —1 .1 l|Gabasch et al.ll200^ . Also, as we de- 
scribe in ii4.5l the p^ measurements can have systematic 
errors at the factor of ^ 2 level which could also account 
for the difference. Obviously resolving the discrepancy 
between these two curves rests as much on an accurate 
determination of the SFR evolution as on that of p*. 

4.4.1. The Stellar Mass Budget 

With the larger number of galaxies afforded by our 
large area we can now split the sample into DRGs and 
non-DRGs to determine their relative contributions to 
the stellar mass budget. Using the mean SEDs for these 
two subsets we calculate A^*/Ly, which is shown in Fig- 
ure The mean A^^/Ly for DRGs is a factor of 1.1 
and 3.6 higher than for the non-DRGs at z = 2.01 and 

"^^ As in R03 both curves have been corrected for the mass loss 
from evolved stars, which asymptotically approaches 3 0% at 13 
Gyr fo r a Salpeter IMF. Changing the IMF to that of IChabrieJ 
(2QP3) would scale both the and SFR(2) measurements down 
by a factor of 1.4 — 1.8. In addition, the mass loss for a Chabrier 
IMF approaches 50% at 13 Gyr and results in a slightly different 
shape of the integral of the SFR(2) curve. 
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2.8 respectively. The higher A^*/L values for DRGs are 
in qualitative agreement with the results from SEP fit- 
ting of individual galaxies ("Forster Schre iber et al.ll2004t 
Ivan Dokkum ct al. 2004; Labbc ct al. 200^! 

We show the corresponding mass densities as the blue 
and red filled pentagons in Figure ^2 The DRGs 
contribute 30% and 64% of the stellar mass density 
at z ~ 2 and 2.8 respectively, compar able to that 
from UV selected samples. As shown in iR^inx et ahl 
( 20031) and more r ecentlv by ' Reddv et all l)2005(l and 
( van Dokkum et al.i .2006) DRGs are almost entirely ab- 
sent in rest-frame UV selected samples like the LEG or 
BM/BX samples. Yet, they make up a comparable frac- 
tion of the mass budget down to similar rest-frame op- 
tical limits, showing that rest-frame UV selection misses 
significant amounts, if not most, of the stellar mass at 
1.6 < z < 3.2. Therefore NIR selection is crucial to ob- 
taining a comprehensive and unbiased view of the high 
redshift Universe. The impor tant contribution of DRGs 
agrees with new results from Ivan Dokkum et all l|200l 
who show that they are the most numerous constituent 
of the population of galaxies with 7V(*> lO^^A^© and 
2 < < 3. It may be that the relatively unobscured star- 
forming UV-selected galaxies are over-represented in our 
luminosity selected sample due to their lower mass-to- 
light ratios. If this is true then the contribution of DRGs 
in a mass-selected sample should therefore be higher. 

Interestingly the mass density of non-DRGs in our 
highest redshift bin is very similar to the HDF-N, which 
is ~ 2 times underdense compared to the total mass 
density from this work and from F03 and R03. This 
may indicate that much of the field-to-field variation at 
high redshift originates in the population of red mas- 
sive galaxies. This was already indicated b y the higher 
clustering amplitude of red galaxies seen bv'Dadd i et al.l 
li2.003) in _thc HDF-S FIRES field and is also found by 
Ivan Dokkum' et all l|2006^) who show that the CDF-S 
is underdense in massive (A^*> 10^^ Mq) galaxies at 
2 < 2 < 3. Clustering measurements of galaxies at high 
redshift as a function of color and over a much larger 
area will directly address this issue (Quadri et al. in 
preparation). 

4.5. Possible Systematic Errors and Biases 

There are various systematic errors which can affect 
our conclusions. First is the limitation that our obser- 
vations only probe out to the rest-frame y-band in our 
highest redshift bin. While stellar mass-to-light ratios 
are much less variable in the rest-frame optical than the 
rest-frame UV they can still vary by an order of mag- 
nitude or more at a given IV, depending on e.g. SFH, 
age, and extinction. Also, the extinction in very dusty 
starbursts can still be quite high in the rest-frame opti- 
cal, implying that the light (and hence mass ) from highly 
extinc ted stars will be missed. The work of lLabbe et al.l 
((200^ and Shaolcv ct al. ( 2005) constrain the possible 
errors by deriving A4i, for objects using rest-frame NIR 
data obtained with the IRAC instrument on the Spitzer 
Space Telescope. For the most part Al* derived from 
these fits agree with those derived from fitting only out 
to the rest-frame optical, although the A^* uncertainties 
are reduced with the longer wavelength baseline. This 
implies that most of the stars that are visible in the rest- 
frame NIR are also visible in the optical, even though 
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Fig. 9. — The evolution of AA-k/hy for luminous galaxies with 
time and redshift, as derived by model fits to the volume averaged 
SEDS. The black pentagons are for the entire lV>L*h''<=sh sample 
and the blue circles and red squares are for non-DRGs and DRGs 
respectively. The open circle represents the local determination 
from the SDSS. 



they may have a higher extinction. If significant stel- 
lar mass exists in environments which are obscured at 
optical and NIR wavelengths then our results would be 
biased. The amount of such heavily obscured stellar mass 
at z < 1 is small, but increa ses rapidly out to higher red- 
shift l)Le Floc'h et al.ll2005ft . If the fraction of extremely 
obscured stellar mass in rest-frame optically luminous 
objects at z > 1 becomes significant then the decline in 
the true mass density will be shallower than what we 
observe. 

Another systematic effect may stem from our relatively 
bright magnitude limits. Even in the HDF-S, which has 
the deepest i^s-band data in existence, we are limited to 
probing high luminosity galaxies and therefore may not 
be measuring a representative sample of the full galaxy 
population. We estimate the amount of light that we 
miss with our IV limit bv using the rest-fram e S-band 
luminosity function of iGiallongo et aP l)2005f) . We con- 
vert their Schechter parameters to the V"-band by cor- 
recting L* using the mean {B — V)rest in each redshift 
bin. At z =0.73, 1.33, 2.01, and 2.8 our luminosity 
limit encompasses approximately 38, 40, 45, and 44% 
of the total V-band light respectively. We also check 
these numbers using the rest-frame y-band luminosity 
function from Marchesini ct al. ( 2006) and find that our 
IV limit encompasses approximately 58 and 70% of the 
light in our two highest redshift bins. This is not to say 
that our relative trends in mass density will be more af- 
fected than other surveys that claim to measure "total" 
mass densities. Some of these surveys, e.g. D03 and 
F03, assume a faint end slope of the luminosity function 
as determined at lower redshift and use it to extrapo- 
late their observed quantities to fainter levels, is then 
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Fig. 10. — The relative evolution of with time and red- 
shift. The symbols are as indicated in the plot. AU quantities 
have been n ormalized to the local stellar mass density measured 
from lCoIe et al. (2001). As described in the text, most other au- 
thors have quoted "total" values determined by extrapolating fits 
to the luminous/massive galaxy population. Our measurements, 
on the other hand, are explicitly restricted to the luminous galaxy 
population, and are complete at 1V>L^'^'''^'^. The dotted errorbars 
on the solid pentagons indicate the minimum and maximum p* 
values among our fields and given an indication of the field-to-field 
variation. The horizontal extent of the solid boxes indicates the 
redshift range over which of our points are derived. The vertical 
extent of the solid boxes indicates the range in values that are 
obtained by allowing the rest-frame luminosity threshold to evolve 
to fainter luminosities at lower redshifts in accordance with our 
derived evolution in A4*/Lv- The mass density averaged over our 
four fields agrees well with most other surveys. The one notable ex- 
ception is the HDF-N (DOS) which has a substantially lower value 
than either the HDF-S or the mean value computed here. The two 
curves are the integrals of the ICole et al.l j2bOM parametric fit to 
the SFR(z) data. The solid magenta line is the fit to the extinction 
corrected data and the dotted blue line line is for the data without 
an extinction correction. As has been noticed by many authors 
an extinction correction to the SFR(2:) measurements is necessary 
to reproduce the p* measurements at almost all redshifts. There 
is a slight systematic difference between the p* estimates and the 
integral of the SFR(z) curve. 



derived by multiplying the resultant luminosity density 
by the mean A^*/L for their directly observed galaxies. 
While formally integrating over all luminosities, the im- 
plicit assumption about the constancy of the faint end of 
the luminosity function is largely uncertain. If the low 
end of the galaxy luminosity function is much steeper 
at high redshifts than at low redshifts this could cause 
the true decline rate in to be less than what we ob- 
serve. In addition, if rest-frame optically faint galaxies 
have very different Ad-^/L values than brighter ones a 
bias in the determinations will occur. The only solu- 
tion to this quandary will be to go significantly deeper 
than the current deepest surveys in the K-hand. Given 
the substantial investment in telescope time required to 
obtain the HDF-S FIRES data (100 hours) only space- 
based observation will be able to push to significantly 



fainter limits. 

In R03 and in Figure H we demonstrated the advan- 
tages of averaging over the galaxy population when de- 
riving stellar mass-to-light ratios. One possible compli- 
cation may result when the different populations being 
averaged have very different extinction properties, such 
as the difference between LBGs and DRGs in our sample. 
We have explored the possible magnitude of such an er- 
ror by fitting various combinations of models with SFHs 
and extinctions that correspond to star- forming DRGs, 
passive DRGs, and LBGs (e.g. POl; Shapley et al. 2001; 
Forster Schreiber et al. 2004; Labbe et al. 2005). Specif- 
ically we have used 3 models consisting of a 1 Gyr old 
CSF model with Ay ^ 2, a 2.9 Gyr old passive pop- 
ulation with no extinction, and a 300 Myr old galaxy 
with Ay ~ 0.6. Depending on the exact contributions 
of the different populations to the mean SED errors of 
up to a factor of 2 in the derived masses can exist. We 
have also checked for the presence of these errors in our 
own data by using the DRG and non-DRG subsamples 
discussed in § 14.4 II On average these two subsamples 
have very different stellar populations and extinctions 
(iForster Schreiber et all 120041: Ivan Dokkum et alJl200l 
iLabbe et al.ll2005|) but the sum of the contributions 
from the two subsamples is within 15% of the value mea- 
sured for the total population, implying that this uncer- 
tainty is not dominant. 

We also have explored the effect of our limited choice 
of SFH on our results. We increased the number of SFHs 

to include SSP, GSF, riooMyr, TsOOMyr, TaOOMyr, T-800Myr, 

■''iGyr, TsGyr: '''eGyr, and a SFH that corres ponds to the 
globally averaged SFR(z) from lCole et"aLl l(200H) . With 
this larger range in SFHs the best-fit p^, values decrease 
by < 15% and the confidence limits in p^, remain un- 
changed. Therefore our choice of SFH is not a large 
source of error. 

Field-to-field variations are an obvious source of error, 
especially since the high luminosity gal axies we are exam- 
ining are likely to b e heavily clustered ijAdelberger et all 
I2004t iDaddi et al]l2003D . We have attempted to miti- 
gate this as much as possible by using data from every 
available field where suitably deep optical and NIR data 
are available. Obviously deep NIR imaging over many 
spatially disjoint fields is crucial to making progress in 
this arena. This will require large investments of time on 
the next generation of wide-field NIR imagers on 8-meter 
class telescopes (e.g. HAWK-I/VLT, MOIRCS/Subaru). 

The IMF is poorly constrained at high redshift. Al- 
though observations of the spectral signatures of massive 
stars in some LBGs indic ate that the h igh mass slope 
must be close to Salpeter ijPettini et al.l l2000) some au- 
thors have argued for the necessity of a top-heavy IMF 
at high redshift to explain abundance ratios in ellipti- 
cal galaxies (e.g. Matteucci 1994 and Nagashima et al. 
2005) and the abundance of submillimete r ga laxies at 
z 2 (Baugh ct al. 200^. As noted in ^ 14.41 changing 
the lower end of the IMF would simply result in a scal- 
ing of all the Pi, and SFR(2;) measurements by the same 
amount. A non-universal IMF, however, either as a func- 
tion of time, environment, or metallicity would result in a 
biased determination of the evolution in and SFR(z). 

In our mod e ling w e assume solar metallicity and a 
iCalzetti et al.l l|200Cl() attenuation law. It is impossible 
to constrain these directly using broad-band photome- 
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try. To determine the sensitivity to these assumptions 
we re-fit our data using Z — O.OO4Z0 models and found 
that Pi, declines by 0.1 — 0.25dex. This is somewhat less 
than the 0.3 — O.Sdex change found by POl and the rea- 
son is not entirely clear. It may be because our mean 
SEDs are significantly redder at all redshifts than the 
very blue starburst SEDs modeled by POl. We explore 
the dependence on the dust atten uation law by fittin g 
using the SMC extinction curve of lPrevot etall l(T98l . 
In this case decreases by < O.ldex. We conclude that 
uncertainties in the metallicity and dust extinction can 
result in up to a 0.2 dex error in our p^, estimates. 

We rely ent irely on BC03 models to interpret the av- 
erage SEDs. Ivan der Wei et all l)2006|) . however, have 
recently pointed out that the BC03 models incorrectly 
predict the evolution in TW^^/L^^ and rest-fra me B — K 
color for early type galaxies at z < 1. The iMarastonI 
()2005(1 models do a better job at the fitting the evolu- 
tion. In the rest-frame optical, the BC03 and the Maras- 
ton models yield identical Aii^/Ly values for all mod- 
els with {U — V)rest > 0.5 but the Maraston models 
yield 7W*/Ly values between 5 and 35% higher for all 
models with < (t7 — V)rest < 0.5 (for solar metallic- 
ity and a Salpeter IMF). Since our (U — V)rest colors 
range from {U — V)rest^ at our highest redshift bin to 
{U — V)rest^ 1 at z = 0, adopting the Maraston mod- 
els would make our mass density evolution shallower by 
5-35%. This is not a dominant source of uncertainty in 
our analysis. 

Our analysis relies heavily on the use of photomet- 
ric redshifts. An incorrect estimate of the redshift will 
lead to errors in the rest-frame colors and luminosities, 
and hence in the stellar masses. Unfortunately the most 
thorough calibrations for Zp^ot estimates have only been 
performed for UV-bright populations for which abun- 
dant spectroscopic redshifts are already in hand. As 
we indicated in the previous section, however, much of 
the stellar mass density at high redshift likely resides 
in UV-faint objects, where Zphot estimates are largely 
untested. Nonetheless, initial progress with optical and 
NIR spectroscopy of DRGs (Kriek et al. 2006; Wuyts et 
al. in preparation) has shown that the photometric red- 
shifts for these objects, while slightly less reliable than 
for optically selected samples, are rarely catastrophically 
wrong. Looking at the whole spectroscopic sample with 
lV>Ly''°'''^ we find that Zphot errors result in systematic 
errors in IV of « 4% and 1% at z < 1.5 and z > 1.5 re- 
spectively. The systematic errors in {U — V)rest are also 
small, !=a 5% and 11%, at z < 1.5 and z > 1.5 respec- 
tively. With the caveat that the spectroscopic samples 
are not representative of the NIR selected galaxy popu- 
lation we conclude that the photometric redshift errors 
should not be a dominant uncertainty in our analysis, al- 
though they are still important for the DRG population. 

A final concern stems from our use of a luminosity as 
opposed to a stellar mass cut to define our galaxy sample. 
Galaxy luminosities can evolve rapidly due to bursts of 
star formation and rapidly changing dust contents. For 
this reason galaxies may pop in and out of a luminosity 
selected sample as they evolve, implying that our sample 
may not be representative of the population as a whole. 
In contrast to the luminosity, A^^, of a galaxy should 
evolve more slowly and a stellar mass selected sample 
will be less susceptible to the above problems. This effect 



may not be significant, since our redshift bins are large 
enough in time to span many typical starburst timescales 
(< a few 100 Myr) and the effect of galaxies entering and 
leaving our luminosity cut may average out. To test the 
dependence of our results on the exact nature of our lu- 
minosity cut we have repeated all measurements using 
a passively evolving luminosity threshold corresponding 
to a maximally old population that has been normal- 
ized to have lV=Ly''°'^'^ at z = 3.2. In the absence of 
dust evolution this cut ensures that all galaxies in our 
high redshift bins will also be present in the lower red- 
shift bins, although galaxies which brightened or stay at 
constant luminosity as a function of time may enter the 
sample at lower redshifts. With this evolving threshold 
the main effect is to increase j"''* and lower 7W^/L as one 
moves to lower redshifts. These changes conspire lower 
p,(z = 2.8)/p*(z = 0) by a factor of 1.9 compared to 
what is shown in Figures [TUl and ITTl and in Tables 13 and 

HI 

We also estimate the change in the inferred evolu- 
tion by evolving our Ly '"^'^ to account for the observed 
change in Aii^/Ly- This translates into a 10.7 times 
fainter luminosity cut at z ~ than at z ^ 2.8. Using 
this fainter threshold, the resultant change in p^, out to 
z ~ 2.8 increases by a factor of 2.1. This would bring our 
measurements further out of agreement with the integral 
of the SFR(z) curve. 

As a related concern, the use of a luminosity limited 
sample also complicates the comparison to other authors 
in § 14.41 Our measured mean A^^/Ly values and our 
j^thrcsh jjj^-^j^ j^i; g^i\ redshifts correspond to a stellar mass 
that is significantly below our mass completeness limit^^. 
This implies that we are likely missing galaxies with high 
Mi,/Lv values. 

In summary, systematic errors in the derived p^, evo- 
lution may be present at the factor of ~ 2 level, with 
the dominant sources of error coming from our lack of 
knowledge about the faint end of the luminosity /mass 
functions and our use of a cut in rest-frame luminosity 
instead of mass. 

5. DISCUSSION 

5.1. The Star Formation Rate Budget 

As discussed in § 14.41 and as shown in Figure E| the 
evolution in pi, roughly a grees with the integral of the 
SFR(z), as calculated bv ICole et all l|200l(l . This SFR 
was determined from rest-frame UV selected samples 
with a large (and uncertain) extinction correction ap- 
plied. That its integral agrees roughly with the direct 
Pi, measurements indicates that the UV-selected samples 
trace most of the SFR density in the Universe, modulo 
the very large extinction corrections. This is in agree- 
ment with recent work of lReddv et al.l l)2005j) and would 
in turn imply that the amount of heavily obscured star 
formation does not dominate the total, as long as the 
stars formed in those heavily obscured environments will 
later be visible in the rest-fra me optical and hence en- 
ter into our p-i, determinations. iReddv et alJ j2005) also 
showed, however, that only < 10% of DRGs would be se- 
lected in traditional UV-selected samples. Coupled with 

■^^ The mass completeness limit that corresponds to our L^'^'^'^'^ 
is derived at each redshift using the A^*/Ly of a maximally old 
single burst with a Salpeter IMF. 
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our measurement of the dominant DRG contribution to 
the mass budget at z > 2 this indicates that UV-selected 
surveys are very incomplete in stellar mass. 

5.2. Constraining the Formation Epoch of Local 
Galaxies 

It is interesting to ask what the observed evolution 
in implies for the formation times of different stellar 
populations in the local Universe. Combining our mea- 
surements with the "total" measurements of other 
authors, the evolution in places the constraint that 
no more than 50% of the stellar mass could have been 
formed at z > 1. Indeed, various authors have found 
that ~ 50% of the l ocal stellar mass in SPSS resides 
in early type ga laxies (' Hogg et a l. 2002: Bell et alJl200.'^ 
iKauffmann et a l. 2003). Age determinations of early 
type galaxies, both from stellar "archeology" (e.g. Trager 
et al. 2000; Thomas et al. 2005), and from fundamen- 
tal plane studies at higher redshifts (e.g. van Dokkum 
et al. 1996; van der Wei et al. 2004) yield formation 
times for the stars in early type galaxies at 1 < z < 3 
in rough agreement with the evolution in p^,. Likewise, 
the remainder of the local stellar mass density must have 
formed at z < 1 compatible with the expected formation 
times for stars in disks of the MW and M31 (e.g. Fergu- 
son & Johnson 2001; Hansen et al. 2002) a nd the inferred 
ages for large disk galaxies at z < 1 ijHammer et al.l 
I200I . 

5.3. Comparison with Theoretical Predictions 

In Figure 1111 we show a comparison of our direct p* 
measurements with the predictions for a set of theoretical 
models. The thick dashed curve shows the evolution in p^ 
derived by integrating an analytic fit to the SFR(z) curve 
as computed from a set of nested hydrodynamic simula- 
tions that include the effects of star formation and feed- 
back (Springel & Hernquist 2003; hereafter SH03). This 
integral takes into account the mass loss from evolved 
stars. The thin solid curve gives the prediction from a 
semi-analytic mock catalog of the Millennium Simulation 
(Springel et al. 2005; Croton et al. 2005; hereafter C05) 
that includes feedback from star formation and from ac- 
tive galactic nuclei (AGN). The total estimates of SH03 
and C05 agree very well with each other, and with the 
local value of p*. 

Because we compute our p* measurements to a fixed IV 
we can compare our data directly to theoretical models 
explicitly including the o bservationa l selection criteria. 
We also include data from lBundv et al. ( 2006) subjected 
to our same IV limit (Bundy private communication). 
Within the substantial ficld-to-field variations (dotted 
errorbars) our p^, est imates and those of [ Bundv et aLI 
(|2OO60 are consistent. iNagamine et al.l 1)20041) claim that 
the discrepancy between total p* estimates and the model 
predictions is caused by observational selection effects. 
According to those authors, the disagreement may be 
caused by the lack of observations at the faint end of 
the galaxy luminosity function coupled with the very 
steep galaxy stellar mas s function in the simula tions. Be- 
cause the simulations of lNagamine et aTl l)2004f) were not 
subjected to observational selection effects, however, it 
is difficult to ascertain the true nature of the disagree- 
ment. We test this explicitly using mock catalogs sub- 
jected to our Ly'''^'^'' limit. The dashed curve is the pre- 



diction of the Millennium simulation. The solid curve 
is derived from a mock galaxy catalog computed from 
the SH03 simulation (Finlator et al. 2006; hereafter 
F06)^^. Both mock catalogs include the effect of extinc- 
ti on on the galaxy luminosities. Contrary to the claim 
of lNagamine et all (^2004) the disagreement between the 
models and the data at high redshift persists even when 
accounting for observational selection in the models. The 
F06 mock catalog prediction has the same general slope 
as the SH03 total p^, estimate and although this model 
comes close to the data at high redshift it fails to match 
the ensemble of data at z < 1.5. The evolution in p^, for 
luminous galaxies in the C05 catalog, however, agrees 
reasonably well with the data at z < 1.5 but produces 
too much mass density in luminous galaxies at z < 1.5. 
In addition, it is worth noting that p^, of luminous galax- 
ies in the Millennium simulation actually decreases from 
z ~ 2 to z ~ 0. 

It is important to note that the overall normalization 
of the simulations is uncertain since the stellar mass de- 
pends intimately on the feedback, which is essentially a 
free parameter of the models, and can accommodate a 
0.3 dex change in the all the curves. The most strik- 
ing disagreement with the models is the large difference 
in the shape of the p^, evolution compared to the obser- 
vations, specifically in that the model evolution is too 
shallow, especially at z < 1.5. 

It is difficult to disentangle the source of the disagree- 
ment between the models and the data. That the model 
curves are too shallow implies that too much mass is 
formed into stars at high redshift, perhaps because of 
an incorrect cooling and feedback prescription. The F06 
and SH03 models do not include AGN feedback which 
may modify the shape. Any such feedback mechanism, 
however, must still be able to reproduce the observed 
number densities of massive galaxies at high redshift (e.g. 
van Dokkum et al. 2006). The slightly decreasing p^, of 
luminous galaxies with cosmic time shown by the C05 
model at z < 2 indicates that the model galaxies are 
evolving rapidly in A^^/L. This cannot be due to the 
mass loss from evolved stars since the C05 simulations 
instantaneously remove 30% of the stellar mass as soon 
as that mass is formed. This is an adequate approxima- 
tion since 25% of the stellar mass is already lost after 
1 Gyr. The most likely explanation for this decline is 
therefore because objects that are above our luminosity 
cut at high redshift fade below it at lower redshifts. The 
AGN feedback prescription adopted by the C05 models is 
the likely cause of this trend as it shuts off star formation 
in massive galaxies, subsequent to which they evolve pas- 
sively to lower redshifts, with the commensurate amount 
of fading. It is puzzling that the models have tradition- 
ally been faulted for not producing the proper number 
of massive galaxies at high redshift but now appear, at 
least for the C05 models, to produce too much mass in 
luminous galaxies at high redshift. 

5.4. Room for Improvement 

There are various ways in which these measurements 
of could be improved, each relating to specific sources 
of possible error: 

Specifically this mock catalog was computed witli the G6 sim- 
ulation. 
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Fig. 11. — A comparison of the evolution of of with the pre- 
dictions of theoretical simulations. The solid pentagons are our 
measurements for galaxies with 1V>L^'^'^'^'', where the dotted er- 
rorbars indicate the minimum and maximum p^. values among our 
fields as an indication of the field-to-field variance. The solid circle 
is our measurement of p* from SDSS for galaxies above the luminos- 
ity limit. The open circle is our estimate of the total p* from SDSS 
and the arrow connecting them indicates our correction to total p* 
for the low red shift point. The o pen hexagon is the total p<: mea- 
surement from Cole et^alj |2001] ). The red empty squares indicate 
the pi, values from Bur Sve^aLr iToog) derived for all galaxies with 
1V>L^''°'''^. The dotted blue curve is the integral of the analytical 
SFR(2) curve from Springel & Hcrnauist (2003), which represents 
the total from that simulation. The solid blue curve is from 
a mock catalog created from the SH03 simulation tl^'inlator et al.l 
120061) and is for simulated galaxies with 1V>L^)?'^''^'^. The results 
of th e Millennium simulation I Springel et al. 2005; Croto n'et al.l 
120051) are plotted as red curves. The dot-dashed red curve cor- 
responds to the total stellar mass density in the simulation and 
the dashed red curve corresponds to the mass density in simulated 
galaxies with lV>Ly^''°'^''. The total p* for all simulations agree 
with the observations at z ~ 0. When compared to the data using 
the same observational cut the p'inlator et al. ( 2006) simulations 
fail to match the ensemble of data at z < 1.4. The Millennium sim- 
ulations, however, fail to match the observational data at z > 1.4. 
We also show the p* values for the non-DRGs and DRGs as filled 
blue circles and red squares respectively, in the two redshift bins 
where DRGs are selected. 



• Field-to-field variations are still a major source of 
uncertainty and can best be addressed by the inclusion 
of many more fields with deep optical and NIR imaging 
data, preferably to K ^ 23. 

• We are still limited in that we are using a rest-frame 
luminosity selected sample. With the deep IRAC data 
that are now becoming available more reliable and well 
constrained mass determinations may be possible for in- 
dividual galaxies. Selection by rest-frame if^-band will 
not entirely solve this problem. I n addition to the m odel 
uncertainties pointed out in M.SllLabbe et alJ l)2005f) and 
IHiiaDlev et al. ( 2005) have demonstrated that there are 
still large (a factor of 6 -13) variations in y \/l^/Lf<- for 
high redshift galaxies and Ivan der Wei et all (2006) have 
shown that A^^^/La' of early type galaxies evolves by a 



factor of 3 over < z < 1. An IRAC-selected sample 
will not be the same as a stellar mass-selected sample. 

• There are few tests of photometric redshifts for galax- 
ies with very red optical-UV colors. Since it is impossi- 
ble to measure spectroscopic redshifts for every object 
it is crucial to test photometric redshifts with extensive 
NIR spectroscopy performed with the new generation of 
multi-object NIR spectrographs. In addition, these ob- 
servations can give us independent measures of the dy- 
namical masses of galaxies and will provide additional 
constraints on the M^, estimates. 

• The uncertainty at the faint end of the luminos- 
ity/mass function is great. Very deep NIR and MIR ob- 
servations will be needed to see if the evolution inferred 
from the rest-frame optically bright population can be 
extended to the population as a whole. 

• Local studies are still inconclusive as to the existence 
of a universal IMF in the local Universe. It is not clear 
at this time how to directly constrain the intermediate 
mass IMF at high redshifts. Dynamical mass estimates 
or gravitational lensing measurements of galaxies in re- 
gions where the baryons are thought to dominate the 
total mass may provide some assistance. 

6. SUMMARY AND CONCLUSIONS 

In this paper we have presented the evolution in the 
properties of the volume averaged stellar population from 
z = to z = 3. We have used deep NIR and optical imag- 
ing over ~ 100 square arcminutes in four disjoint fields to 
mitigate field-to-field variations. The deep NIR data al- 
low us to select and analyze galaxies at rest-frame optical 
wavelengths, where dust extinction and stellar mass-to- 
light ratio variations are far less than in the rest-frame 
UV. We derived photometric redshifts for all objects lack- 
ing a spectroscopic redshift and measure rest-frame UV 
through optical luminosities of all galaxies at z < 3.2. 
We create a rest-frame F-band luminosity selected sam- 
ple with 1V> 3 X 10^° h^Q^ Lq above which we are com- 
plete at all redshifts. 

From these individual measurements we construct vol- 
ume averaged properties such as the luminosity density, 
and color, e.g. {U - B)rest and {B - V)rest- The 
redshift evolution in luminosity density and mean color 
is very different. As found in R03 for a much smaller 
sample, jy^^ is roughly constant within a factor of « 2.2 
from z ~ 3 to z = 0. This modest evol ution is broadly 
consis tent with the results of D03 and iGiallongo et alJ 
|)200!t^ . The evolution in the rest-frame colors, however, 
is dramatic and the colors become systematically red- 
der with decreasing redshift. The colors of our SDSS 
comparison sample at z = 0.1 are {B — V)rest~ 0.72 
and {U — B)rest= 0.07. In contrast the average color 
of luminous galaxies at z ~ 3 is (-B — V)rest= 0.40 and 

{U~B)rest^ -0.38. 

By averaging over the population as a whole we are 
also averaging over the individual SFHs of the galaxies 
and the resultant volume averaged SFH will be smoother 
than that of the individual galaxies. We expect the mean 
colors to be more easily modeled than those of individual 
galaxies, which presumably have bursty SFHs. We find 
that the mean colors at all redshifts lie close to the locus 
of morphologically normal local galaxies in {B~V)rest vs. 
([/ — B)rest space but compared to the local sample are 
slightly bluer in {U — B)rest at a fixed [B — V)rest- Some 



Rudnick et al. 



19 



of this difference could result from photometric redshift 
errors but could also be indicative of a residual signature 
of bursts in the mean colors. 

We derive the luminosity weighted volume averaged 
UV-optical rest-frame SED for all galaxies with 1V> 
3 X 10^^ hjQ Lq. At all wavelengths the mean SED be- 
comes gradually more red with decreasing redshift and 
at all redshifts is within the range of SEDs spanned by 
normal galaxies in the local Universe. Even at our high- 
est redshift bin the volume average SED d oes not look 
like that of a starburst from Kinnev et all fl996 ) with 
little extinction but rather, is redder than an irregular 
template. We also compute the mean SED from DRGs 
chosen to have (Js - Ks)> 2.3. This mean DRG SED 
redder than average in the rest-frame optical and in the 
rest-frame UV. 

With simple models we transform the mean SED into 
an evolution in the global yW*/Ly. Because the volume 
averaged SFH should be smoother than that of the in- 
dividual galaxies it is therefore more appropriate to use 
smooth models to fit the mean SED than to fit the SEDs 
of individual galaxies. Nonetheless, the colors cannot be 
fit well at all redshifts with a single model having an 
exponential SFH and a constant dust obscuration and 
metallicity. We can, however, fit the average SED at 
each redshift reasonably well with some combination of 
the parameters, as long as they can change as a function 
of redshift. We use these model fits to derive luminosity 
weighted stellar mass-to-light ratios for all galaxies with 
1V> 3 X 10^" h^Q Lq. Mi,/Lv also declines smoothly 
from z = 0.1 to z ~ 3 with A^^/Ly= 4.9 and 0.5 at 
z = 0.1 and z ~ 3 respectively. There are strong varia- 
tions in Aii,/Lv across the population. The DRGs have 
an Ai^/Lv a factor of « 1.1 and 3.6 higher than the non- 
DRGs in the z = 2.01 and z = 2.8 redshift bins. This 
is consistent with the r e sults of Jjabbc ct al. ( 2005) and 
iForster Schreiber et all l)2004() who modeled the SEDs of 
individual DRGs. 

Multiplying the jy^^ and mean A^*/Ly measurements 
we derive the evolution in out to z '--^ 3. in luminous 
galaxies declines by a factor of 4.8 out to z = 3. At 
z - 2 and 2.8, the DRGs contribute 30% and 64% of 
Pi, showing that UV selection techniques miss much of 
the stellar mass at 1.6 < z < 3 in rest-frame optically 
luminous galaxies. 

We also compare our p^, measurements to a set of re- 
cent model predictions subjected to our same rest-frame 
luminosity cut and find a large disagreement between 
the models themselves and between the models and ob- 
servations. The different models can match the observed 
evolution in p^ over a limited redshift range, but both 
models predict an evolution that is too shallow to match 
the observations. This therefore implies that the masses 
of luminous galaxies in the models, or their abundance, 
are not consistent with the observations. 

To better understand how different types of galaxies 
contribute to the stellar mass budget at high redshifts it 
will be necessary to reduce our reliance on photometric 
redshifts by performing NIR spectroscopy on a represen- 
tative sample of DRGs. This will confirm the redshifts 
of the optically faint population that contributes over 
50% of the stellar mass at high redshift. Re-performing 
this experiment with IRAC data will reduce the uncer- 
tainty due to dust obscuration and obtaining very deep 



NIR/MIR imaging will allow us to constrain the low mass 
end of the galaxy stellar mass function. 
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TABLE 1 

Ultra-violet Rest-frame Luminosity Densities 



z 


Fiold(s) 


Ntot 


^use 


32200 


4?2200 




4?2700 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


0.73 


ALL 


116 


112 


-11.84 


0.34 


-12.61 


0.14 


1.33 


ALL 


154 


127 


-12.35 


0.18 


-12.68 


0.14 


2.01 


ALL 


164 


126 


-11.95 


0.14 


-12.38 


0.13 


2.80 


ALL 


95 


82 


-13.67 


0.15 


-13.84 


0.14 


0.73 


HDF-N 


11 


11 


-12.06 


0.42 


-12.79 


0.41 


1.33 


HDF-N 


9 


9 


-13.15 


0.55 


-13.46 


0.50 


2.01 


HDF-N 


13 


11 


-12.91 


0.38 


-13.18 


0.37 


2.80 


HDF-N 


16 


16 


-13.82 


0.30 


-14.02 


0.30 


0.73 


HDF-S 


4 


4 


-11.55 


0.58 


-12.12 


0.60 


1.33 


HDF-S 


13 


11 


-12.78 


0.33 


-13.22 


0.33 


2.01 


HDF-S 


7 


6 


-12.55 


0.51 


-12.84 


0.50 


2.80 


HDF-S 


25 


24 


-14.60 


0.25 


-14.73 


0.24 


1.33 


MS1054-03 


36 


34 


-12.64 


0.24 


-12.90 


0.22 


2.01 


MS1054-03 


50 


46 


-12.35 


0.25 


-12.77 


0.24 


2.80 


MS1054-03 


54 


42 


-13.27 


0.22 


-13.46 


0.20 


0.73 


CDF-S 


101 


97 






-12.62 


0.15 


1.33 


CDF-S 


96 


73 


-12.07 


0.29 


-12.44 


0.22 


2.01 


CDF-S 


94 


63 


-11.54 


0.20 


-12.04 


0.18 



Note. — All luminosity densities are given in /170 AB magnitudes / 
Mpc^. (3) is the total number of galaxies in each bin with 1V>L^'^°'''^. 
In (4) we give the number of galaxies above our lum inos ity threshold 
that also have Sz/{1 + z) < 0.16. As described in i|3.2l and R03 we 
calculate the j'^'^'^* values using the galaxies in (4), and correct for those 
galaxies that have been excluded due to their large photometric redshift 
uncertainties. The z = 0.73 entry for MS1054-03 has been omitted 
because of the presence of a rich 2 = 0.83 cluster in that field. The 
z = 2.80 entry for the CDF-S has been omitted because of a high 
incompleteness at these redshifts caused by the moderate depth. The 
empty entry in columns (5) and (6) for the CDF-S is because Xrest = 
22OOA was blueward of the bluest observed filter. 



TABLE 2 

Optical Rest-frame Luminosity Densities 



z 


Field (s) 


ju 




dju 








B 


jv 




djv 


(1) 


(2) 


(3) 




(4) 


(5) 




(6) 


(7) 




(8) 


0.73 


ALL 


-13, 


,86 


0.12 


-14, 


,82 


0, 


,12 


-15 


,40 


0.12 


1.33 


ALL 


-13, 


,73 


0.09 


-14, 


,52 


0, 


,09 


-14, 


.98 


0.09 


2.01 


ALL 


-13, 


,53 


0.09 


-14, 


,31 


0, 


,09 


-14, 


,78 


0.08 


2.80 


ALL 


-14, 


,55 


0.11 


-15, 


,05 


0, 


,12 


-15, 


,35 


0.12 


0.73 


HDF-N 


-14, 


12 


0.40 


-15, 


,15 


0, 


,38 


-15, 


,73 


0.37 


1.33 


HDF-N 


-14, 


,26 


0.41 


-15, 


,03 


0, 


,39 


-15, 


,52 


0.41 


2.01 


HDF-N 


-14, 


,01 


0.28 


-14, 


,59 


0, 


,28 


-14, 


,93 


0.30 


2.80 


HDF-N 


-14, 


,48 


0.31 


-14, 


,90 


0, 


,29 


-15 


.20 


0.28 


0.73 


HDF-S 


-13, 


,21 


0.59 


-14, 


,12 


0, 


,61 


-14, 


,64 


0.62 


1.33 


HDF-S 


-14, 


,36 


0.26 


-15, 


,17 


0, 


,27 


-15, 


,62 


0.26 


2.01 


HDF-S 


-13, 


,63 


0.40 


-14, 


,20 


0, 


,38 


-14, 


,46 


0.38 


2.80 


HDF-S 


-15, 


,22 


0.21 


-15, 


,70 


0, 


,19 


-15, 


,93 


0.19 


1.33 


MS1054-03 


-13, 


,72 


0.18 


-14, 


,43 


0, 


,17 


-14, 


,85 


0.17 


2.01 


MS1054-03 


-13, 


,76 


0.18 


-14, 


,55 


0, 


,17 


-15, 


,00 


0.17 


2.80 


MS1054-03 


-14, 


,35 


0.15 


-14, 


,89 


0, 


,16 


-15, 


,21 


0.17 


0.73 


CDF-S 


-13, 


,87 


0.13 


-14, 


,83 


0, 


,13 


-15, 


,40 


0.12 


1.33 


CDF-S 


-13, 


,62 


0.12 


-14, 


,44 


0, 


,11 


-14, 


,91 


0.11 


2.01 


CDF-S 


-13, 


,37 


0.11 


-14, 


,18 


0, 


,11 


-14, 


,69 


0.11 



Note. — All luminosity densities are given in ^70 AB magnitudes 
/ Mpc^ with the following conversion to Vega magnitudes: Uvega = 
Uab-0.79, B^ega = Bas+0.102, Vvega = Vas-0.008. The z = 0.73 
entry for MS1054-03 has been omitted because of the presence of a 
rich z = 0.83 cluster in that field. The z = 2.80 entry for the CDF-S 
has been omitted because of a high incompleteness at these redshifts 
caused by the moderate depth. 
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TABLE 3 

Stellar Mass Densities 



V-*-/ 


T thresh 
(2) 


log(p*) 


f41 
V^} 


(5) 


(61 


f7) 


0.1 





8.59 


8.55 


8.63 






0.1 


3 X IQio 


8.27 


8.23 


8.31 






0.73 


3 X IQi" 


7.97 


7.83 


8.03 


7.63 


8.52 


1.33 


3 X IQio 


7.81 


7.74 


7.87 


7.63 


8.12 


2.01 


3 X 10^° 


7.73 


7.70 


7.82 


7.08 


7.82 


2.80 


3 X IQio 


7.59 


7.47 


7.67 


7.38 


7.73 


0.1 


0.3 X lOio ^ 


8.59 


8.55 


8.63 






0.73 


1.3 X 10^° ^ 


8.05 


8.02 


8.12 


7.86 


8.52 


1.33 


1.3 X 10^° * 


7.87 


7.83 


7.94 


7.74 


8.22 


2.01 


1.3 X 10i° * 


7.76 


7.70 


7.82 


7.31 


7.89 


2.80 


3 X 10i° ^ 


7.59 


7.48 


7.65 


7.38 


7.73 



Note. — (2) is given in units of h^Q^ Ly,©. (3)-(7) are given in units of hyo M.Q Mpc ^ 
and are for all galaxies with 1V>L*''*'=''. The first 2 rows indicate our measurements for the 
SDSS. (4) and (5) are the formal confidence limits on p*. (6) and (7) give the range in 
allowed by the rms field-to-field variations among our different fields. 

''These IV limits are computed using the observed variation in Jvl^,/hv for galaxies with 
1V> 3 X lOiOh-2 Lv,©. 



TABLE 4 
Mean Mass-to-light ratios 



z 


Field(s) 


(M,/Lv) 


{Mi,/Lv)iow 


{Mi, /Lv) high 


(1) 


(2) 


(3) 


(4) 


(5) 


0.73 


ALL 


1.06 


0.77 


1.18 


1.33 


ALL 


1.07 


0.88 


1.18 


2.01 


ALL 


1.08 


0.97 


1.19 


2.80 


ALL 


0.46 


0.39 


0.50 


0.73 


HDF-N 


1.30 


1.06 


1.44 


1.33 


HDF-N 


0.97 


0.65 


1.18 


2.01 


HDF-N 


0.58 


0.40 


0.76 


2.80 


HDF-N 


0.33 


0.15 


0.39 


0.73 


HDF-S 


0.96 


0.85 


1.36 


1.33 


HDF-S 


1.19 


0.78 


1.26 


2.01 


HDF-S 


0.32 


0.28 


0.62 


2.80 


HDF-S 


0.37 


0.22 


0.43 


1.33 


MS1054-03 


0.80 


0.71 


0.97 


2.01 


MS1054-03 


1.08 


0.78 


1.19 


2.80 


MS1054^03 


0.50 


0.37 


0.62 


0.73 


CDF-S 


1.06 


0.96 


1.06 


1.33 


CDF-S 


1.08 


0.78 


1.18 


2.01 


CDF-S 


1.26 


1.03 


1.39 



Note. — Columns (3)- (5) arc computed for all galaxies with 
lV>L*h''<=* and arc in units of [Mq/Lq]. The z = 0.73 entry for 
MS1054-03 has been omitted because of the presence of a rich 
2 = 0.83 cluster in that field. The z = 2.80 entry for the CDF- 
S has been omitted because of a high incompleteness at these 
redshifts caused by the moderate depth. 
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TABLE 5 
Relative Stellar Mass Densities 



2 
(1) 


Fiold(s) 
(2) 


(3) 


Ap*,!otu 

(4) 


(5) 


/^pcorr 

-k.iOW 

(6) 


^pThLh 

(7) 


0.73 


ALL 


U 




0.28 


0.44 


U 


.AL 


0.44 


1.33 


ALL 


U 


.A i 


0.23 


0.30 


U 


1 Q 


0.30 


2.01 


ALL 


U 




0.21 


0.27 


u. 


1 n 
. iu 


0.27 


2.80 


ALL 


U 


1 P. 


0.12 


0.19 


U. 


.UO 


0.19 


0.73 


HDF-N 


0, 


.64 


0.38 


0.83 


0, 


.18 


0.83 


1.33 


HDF-N 


0, 


.39 


0.19 


0.57 


0, 


.09 


0.57 


2.01 


HDF-N 


0, 


.14 


0.08 


0.22 


0, 


.04 


0.22 


2.80 


HDF-N 


0, 


.10 


0.05 


0.13 


0, 


.02 


0.13 


0.73 


HDF-S 


0, 


.17 


0.07 


0.32 


0, 


.04 


0.32 


1.33 


HDF-S 


0, 


.53 


0.33 


0.69 


0, 


.16 


0.69 


2.01 


HDF-S 


0, 


.05 


0.04 


0.10 


0, 


.02 


0.10 


2.80 


HDF-S 


0, 


.22 


0.13 


0.28 


0, 


.06 


0.28 


1.33 


MS1054-03 


0, 


.18 


0.14 


0.22 


0, 


.07 


0.22 


2.01 


MS1054-03 


0, 


.27 


0.19 


0.33 


0, 


.09 


0.33 


2.80 


MS1054-03 


0, 


.15 


0.10 


0.21 


0, 


.05 


0.21 


0.73 


CDF-S 


0, 


.39 


0.34 


0.44 


0, 


.16 


0.44 


1.33 


CDF-S 


0, 


.25 


0.19 


0.29 


0, 


.09 


0.29 


2.01 


CDF-S 


0, 


.24 


0.20 


0.29 


0, 


.10 


0.29 



Note. — The relative mass densities in (3) are given in terms of the 
fractional decrease from z = 0, i.e. Ap* = pt/p*(2 = 0), as computed 
by correcting the valu es to total using t he SDSS data and using the 
2 = measurement of ICole et al.l 1200 U) . The confidence intervals in 
Api given in columns (4) and (5) come from the formal uncertainties 
on this evolution for galaxies with 1V>L^'''''''^. The confidence intervals 
in columns (6) and (7) include the extra uncertainty that comes from 
allowing L^'""'' to evolve with redshift according to the evolution in 
A^*/Lv- The 2 = 0.73 entry for MS1054-03 has been omitted because 
of the presence of a rich z = 0.83 cluster in that field. The z = 2.80 
entry for the CDF-S has been omitted because of a high incompleteness 
at these redshifts caused by the moderate depth. 
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APPENDIX 

THE EFFECT OF WEAK LENSING ON THE LUMINOSITY DENSITY 

Weak lensing of background galaxies by a foreground mass distribution affects the luminosity density of the back- 
ground sources by changing their apparent magnitudes, surface densities, and the volume enclosed in the observed 
solid angle. To the first order all of these effects are directly dependent on the magnification factor ji. Given a mass 
map for the foreground cluster and the redshift of each source it is possible to compute the magnifications of each 
source fii. 

First we consider the effect of weak lensing on the total integrated luminosity density in a given redshift slice whose 
enclosed co-moving volume per unit solid angle is Vq, in units of [Mpc'^/ster]. The "true" luminosity density, i.e. in 
the absence of lensing, is jo- This can be written as: 



^gaL 

no Lo^i 
Vq ^ n„ 

4 = 



Mpc3 



i-n '"gal 

where uq is the surface density of sources in units of [ster~^]. ng can be written as 

no = NoVo [ster-i], 

where A'o is the intrinsic volume number density in units of [^/Mpc^]. We can now rewrite jo as 



30 



2 — 



[Lo/Mpc^ 



(Al) 



(A2) 



(A3) 



Under the influence of lensing the rest-frame luminosities (which is what we derive from observed magnitudes) 
become 



(A4) 



while the intrinsic volume number density naturally remains unaffected. The luminosity density subject to lensing 
ji is therefore 



ji = No 



"soi j- 



2—0 ^ 



(A5) 



Since the sole difference between jg and jo is the magnification of the individual sources, we compute jo in MS1054-03 
by using the demagnified galaxy luminosities, Loa = Lg^ij [li. The /ii are determined from the weak lensing map of 
iHockstra et al. (2000). 

The second effect stems from the use of a fixed luminosity limit in the presence of magnification. For example, 
sources whose intrinsic luminosities would exclude them from the sample can enter the sample after being magnified. 
At all but the highest redshifts in the MS1054-03 field our data are complete significantly below the luminosity limit. 
Therefore to correct for this effect in our analysis we select galaxies brighter than L^''°''^ based on their demagnified 
luminosities. 



